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Preface 



The motion of classical mechanical systems is determined by Hamilton’s dif- 
ferential equations: 

(x(t) = d y H(x(t),y(t)) 

\j/(t) = -d x H(x(t),y(t)) 

For instance, if we consider the motion of n particles in a potential field, the 
Hamiltonian function 

1 n 

H = V ( x l !••••> x n ) 

Z 2 = 1 

is the sum of kinetic and potential energy; this is just another formulation of 
Newton’s Second Law. 

A distinguished class of Hamiltonians on a cotangent bundle T*X con- 
sists of those satisfying the Legendre condition. These Hamiltonians are ob- 
tained from Lagrangian systems on the configuration space X , with coordi- 
nates (x,x) = (space, velocity ), by introducing the new coordinates (x, y) = 
(space, momentum) on its phase space T*X. Analytically, the Legendre con- 
dition corresponds to the convexity of H with respect to the fiber variable y. 
The Hamiltonian gives the energy value along a solution (which is preserved 
for time-independent systems) whereas the Lagrangian describes the action. 
Hamilton’s equations are equivalent to the Euler-Lagrange equations for the 
Lagrangian: 

^-d±L(x(t),x(t)) = d x L(x(t),x(t)). 
at 

These equations express the variational character of solutions of the La- 
grangian system. A curve x : [to? ^l] — » K" is a Euler-Lagrange trajectory 
if, and only if, the first variation of the action integral, with end points held 
fixed, vanishes: 

r* 1 x(ti) 

S / L(x(t), x(t)) dt =0. 

Jt 0 -*(*o) 
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In other words, solutions extremize the action with fixed end points on each 
finite time interval. 

This is not quite what one usually remembers from school 1 , namely that 
solutions should minimize the action. The crucial point here is that the min- 
imizing property holds only for short times. For instance, when looking at 
geodesics on the round sphere, the movement along a great circle ceases to be 
the shortest connection as soon as one comes across the antipodal point. 

However, under certain circumstances there may well be action minimizing 
trajectories. The investigation of these minimal objects is one of the central 
topics of the present work. In fact, they do not always exist as genuine solu- 
tions, but they do so as invariant measures. This is the outcome of a theory by 
Mather and Mane which generalizes Aubry-Mather theory from one to more 
degrees of freedom. In particular, there exist action minimizing measures with 
any prescribed “asymptotic direction” (described by a homological rotation 
vector) . Associating to each such rotation vector the action of a minimal mea- 
sure, we obtain the minimal action functional 

ex : F/i(A, M) — * M. 

By construction, the minimal action does not describe the full dynamics but 
concentrates on a very special part of it. The fundamental question is how 
much information about the original system is contained in the minimal ac- 
tion? 

The first two chapters of this book provide the necessary background on 
Aubry-Mather and Matlrer-Mane theories. In the following chapters, we in- 
vestigate the minimal action in four different settings: 

1. convex billiards 

2. fixed points and invariant tori 

3. Hofer’s geometry 

4. symplectic geometry. 

We will see that the minimal action plays an important role in all four situa- 
tions, underlining the significance of that particular variational principle. 

1. Convex billards. Can one hear the shape of a drum? This was Kac’ pointed 
formulation of the inverse spectral problem: is a manifold uniquely determined 
by its Laplace spectrum? We do know now that this is not true in full gen- 
erality; for the class of smooth convex domains in the plane, however, this 
problem is still open. 

We ask a somewhat weaker question for the length spectrum (i.e., the set 
of lengths of closed geodesics) rather than the Laplace spectrum, which is 
closely related to the previous one: how much geometry of a convex domain 
is determined by its length spectrum? The crucial observation is that one can 
consider this geometric problem from a more dynamical viewpoint. Namely, 

1 depending on the school, of course. . . 




Preface VII 



following a geodesic inside a convex domain that gets reflected at the bound- 
ary, is equivalent to iterating the so-called billiard ball map. The latter is a 
monotone twist map for which the minimal action is defined. 

The main results from Chapter 3 can be summarized as follows. 

Theorem 1 . For planar convex domains, the minimal action is invariant un- 
der continuous deformations of the domain that preserve the length spectrum. 

In particular, every geometric quantity that can be written in terms of the 
minimal action is automatically a length spectrum invariant. 

In fact, the minimal action is a complete invariant and puts all previously 
known ones (e.g., those constructed in [2, 19, 63, 87]) into a common frame- 
work. 

2. Fixed points and invariant tori. We consider a symplectic diffeomorphism 
in a neighbourhood of an elliptic fixed point in R 2 . If the fixed point is of 
“general” type, the symplectic character of the map makes it possible (under 
certain restrictions) to find new symplectic coordinates in which the map 
takes a particularly simple form, the so-called Birklroff normal form. The 
coefficients of this normal form, called Birklroff invariants, are symplectically 
invariant. 

The Birklroff normal form describes an asymptotic approximation, in the 
sense that it coincides with the original map only up to a term that vanishes 
asymptotically when one approaches the fixed point. In general, it does not 
give any information about the dynamics away from the fixed point. 

The main result in this context introduces the minimal action as a sym- 
plectically invariant function that contains the Birklroff normal form, but also 
reflects part of the dynamics near the fixed point. 

Theorem 2. Associated to an area-preserving map near a general elliptic 
fixed point there is the minimal action a, which is symplectically invariant. 

It is a local invariant, i.e., it contains information about the dynamics 
near the fixed point. Moreover, the Taylor coefficients of the convex conjugate 
a* are the Birkhoff invariants. 

Area-preserving maps near a fixed point occur as Poincare maps of closed 
characteristics of three-dimensional contact flows. A particular example is 
given by the geodesic flow on a two-dimensional Riemannian manifold. In 
this case, the minimal action is determined by the length spectrum of the 
surface, and we obtain the following result. 

Theorem 3. Associated to a general elliptic closed geodesic on a two-dimen- 
sional Riemannian manifold there is the germ of the minimal action, which is 
a length spectrum invariant under continuous deformations of the Riemannian 
metric. 

The minimal action carries information about the geodesic flow near the 
closed geodesic; in particular, it determines its C° -integrability . 
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In higher dimensions, we consider a symplectic diffeomorplrism <f> in a 
neighbourhood of an invariant torus A. If we assume that the dynamics on A 
satisfy a certain non-resonance condition, one can transform <j> into Birklroff 
normal form again. If this normal form is positive definite the map </> deter- 
mines the germ of the minimal action a, and we will show again that the 
minimal action contains the Birklroff invariants as Taylor coefficients of a*. 

3. Hofer’s geometry. Whereas the first three settings had many features in 
common, the viewpoint here is quite different. Instead of looking at a single 
Hamiltonian system, we investigate all Hamiltonian systems on a symplectic 
manifold (M, ui) at once, collected in the Hamiltonian diffeomorplrism group 
Ham(M, u>). It is one of the cornerstones of symplectic topology that this group 
carries a bi-invariant Finsler metric d, usually called Hofer metric, which is 
constructed as follows. 

Think of Ham(M, lo) as infinite-dimensional Lie group whose Lie algebra 
consists of all smooth, compactly supported functions H : M — ■> R. with mean 
value zero. Introduce any norm || • || on those functions that is invariant under 
the adjoint action H i— > Hoip- 1 . Then the Hofer distance of a diffeomorplrism 
(f) from the identity is defined as the infrmum of the lengths of all paths in 
Ham(M, u>) that connect 4> to the identity: 

d(id, <f>) = inf { [ \\H t \\ dt \ gA H = (j)}. 

Jo 

The problem is to choose the norm || • ||. The Hamiltonian system is determined 
by the first derivatives of H , but ||di?j|co, for instance, is not invariant under 
the adjoint action. It turns out that the oscillation norm 



|| • || = osc := max — min 



is the right choice although it seems to have nothing to do with the dynamics. 
Loosely speaking, the Hofer metric generates a C _1 -topology and measures 
how much energy is needed to generate a given map. 

The resulting geometry is far from being understood completely. This is 
due to the fact that, despite its simple definition, the Hofer distance is very 
hard to compute. After all, one has to take all Hamiltonians into account 
that generate the same time-l-map. A fundamental question concerns the re- 
lation between the Hofer geometry and dynamical properties of a Hamiltonian 
diffeomorplrism: does the dynamical behaviour influence the Hofer geometry 
and, vice versa, can one infer knowledge about the dynamics from Hofer’s 
geometry? Only little is known in this direction. 

In Chap. 5, we take up this question for Hamiltonians on the cotan- 
gent bundle T*T” satisfying a Legendre condition. This leads to convex La- 
grangians on T T" for which the minimal action a is defined. On the other 
hand, the Hamiltonians under consideration are unbounded and do not fit 
into the framework of Hofer’s metric. Therefore, we have to restrict them to 
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a compact part of T* T", e.g., to the unit ball cotangent bundle B* T™, but in 
such a way that we stay in the range of Mather’s theory. 

Let a denote the minimal action associated to a convex Hamiltonian diffeo- 
morphism on B*T n . Our main result in this context shows that the oscillation 
of a*, which is nothing but a(0), is a lower bound for the Hofer distance. This 
establishes a link between Hofer’s geometry of convex Hamiltonian mappings 
and their dynamical behaviour. 

Theorem 4 . If <f> G Ham(B* T") is generated by a convex Hamiltonian then 

d(id,(/)) > osc a* = a(0). 

f. Symplectic geometry. Consider the cotangent bundle T* T" with its canon- 
ical symplectic form u > 0 = dA. Here, A is the Liouville 1-form which is ydx in 
local coordinates ( x , y). Suppose H : T* T" — » R. is a convex Hamiltonian. Be- 
cause H is time-independent the energy is preserved under the corresponding 
flow, i.e. , all trajectories lie on (hberwise) convex (2 n — l)-dimensional hyper- 
surfaces £ = {H = const.}. Of particular importance in classical mechanics 
are so-called KAM-tori. i.e., invariant tori carrying quasiperiodic motion. 
These are graphs over the base manifold T", with the additional property 
that the symplectic form uj 0 vanishes on them; submanifolds with the latter 
property are called Lagrangian submanifolds. 

We want to study symplectic properties of Lagrangian submanifolds on 
convex hypersurfaces. To do so, we observe that a Lagrangian submanifold 
possesses a Liouville class a a, induced by the pull-back of the Liouville form 
A to A. The Liouville class is invariant under Hamiltonian diffeomorphisms, 
i.e., it belongs to the realm of symplectic geometry. On the other hand, be- 
ing a graph is certainly not a symplectic property. Our starting question in 
this context is as follows: is it possible to move a Lagrangian submanifold A 
on some convex hypersurface A by a Hamiltonian diffeomorphism inside the 
domain Us bounded by £1 

In a first part, we will see that, under certain conditions on the dynamics 
on A, it is impossible to move A at all; we call this phenomenon boundary 
rigidity. In fact, the Liouville class a a already determines A uniquely. 

Theorem 5. Let A be a Lagrangian submanifold with conservative dynamics 
that is contained in a convex hypersurface £, and let K be another Lagrangian 
submanifold inside Us- Then 



a a = ok A = K. 

What can happen if boundary rigidity fails? Surprisingly, even if it is pos- 
sible to push A partly inside the domain Us, it cannot be done completely. 
Certain pieces of A have to stay put, and we call them non-removable inter- 
sections. In the case where £ is some distinguished “critical” level set, these 
non-removable intersections always contain an invariant subset with specific 
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dynamical behaviour; this subset is the so-called Aubry set from Mather 
Mane theory. This result reveals a hidden link between aspects of symplectic 
geometry and Mather-Mane theory in modern dynamical systems. 

Finally, we come back to the somewhat annoying fact that the property 
of being a Lagrangian section is not preserved under Hamiltonian diffeomor- 
phisms. For this, we consider 

Theorem 6. Let U be a (fiberwise) convex subset U of T* T". Then every 
cohomology class that can be represented as the Liouville class of some La- 
grangian submanifold in U, can actually be represented by a Lagrangian sec- 
tion contained in U . 

So, from this rather vague point of view at least, Lagrangian sections actually 
do belong to symplectic geometry. 

Furthermore, the above result allows symplectic descriptions of seemingly 
non-symplectic objects: the stable norm from geometric measure theory, and 
also our favourite, the minimal action. 

Theorem 7. The stable norm of a Riemannian metric g on T™ , and the min- 
imal action of a convex Lagrangian L : T T™ — > R, both admit a symplectically 
invariant description. 

This closes the circle for our investigation of the Principle of Least Action 
in geometry and dynamics. 
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Aubry— Mather theory 



The Principle of Least Action states that, for sufficiently short times, tra- 
jectories of a Lagrangian system minimize the action amongst all paths in 
configuration space with the same end points. If the time interval becomes 
larger, however, the Euler-Lagrange equations describe just critical points of 
the action functional; they may well be saddle points. 

In the eighties, Aubry [5] and Mather [64] discovered independently that 
monotone twist maps on an annulus possess orbits of any given rotation num- 
ber which minimize the (discrete) action with fixed end points on all time 
intervals. Roughly speaking, the rotation number of a geodesic describes the 
direction in which the geodesic, lifted to the universal cover, travels. Those 
minimal orbits turned out to be of crucial importance for a deeper under- 
standing of the complicated orbit structure of monotone twist mappings. 

Later, Mather [69] developed a similar theory for Lagrangian systems in 
higher dimensions. There was, however, an old example by Hedlund [41] of 
a Riemannian metric on T 3 , having only three directions for which minimal 
geodesics existed. Therefore, Mather’s generalization deals with minimal in- 
variant measures instead of minimal orbits. 

A different approach was suggested by Mane [62] who introduced a certain 
critical energy value at which the dynamics of a Lagrangian systems change. 
It turned out that this approach essentially contains Mather’s theory, but in 
a more both geometrical and dynamical setting. 

We will deal with these generalizations of Aubry-Matlrer theory to higher 
dimensions in Chap. 2. 



1.1 Monotone twist mappings 

Let 

S 1 x (a, b) C S 1 xR = rs' 

be a plane annulus with S 1 = R/Z, where we allow the cases a = — oo or 
b = +oo (or both). Given a diffeomorphism (j> of S 1 x (a, b) we consider a lift (f> 
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1 Aubry-Mather theory 



of (j) to the universal cover R x (a, b) of S 1 x (a, b ) with coordinates x, y. Since 
0 is a diffeomorphism, so is </>, and we have <p{x + 1 , y) = cf>(x,y) + (1,0). In 
this section, we will always work with (fixed) lifts for which we drop the tilde 
again and keep the notation </>. 

In the case when a or & is finite we assume that f> extends continuously to 
R x [a, b } by rotations by some fixed angles: 

<fi(x, a) = (x + u>-, a) and f>(x, b) = (x + w+, b). (1-1) 



The numbers u>± are unique after we have fixed the lift. For simplicity, we set 
u>± = ±oo if a = —oo or 6 = oo. 

Definition 1 . 1 . 1 . A monotone twist map is a C 1 -diffeomorphism 



(f> : R x (a, b) — » R x (a, b) 

(x 0 ,yo) (xi,yi) 



satisfying f>{x o + 1, yo) = <f>{ xq, yo) + (1, 0) as well as the following conditions: 

1. 4> preserves orientation and the boundaries of R x ( a,b ), in the sense that 
Vi(xo,yo) a,b as y 0 -» a, b; 

2. if a or b is finite <f> extends to the boundary by a rotation, i.e., it satisfies 

(i-i); 

3. 4> satisfies a monotone twist condition 




(1.2) 



4- 4> is exact symplectic; in other words, there is a C 2 -function h, called a 
generating function for <p, such that 



yidxi - y 0 dx 0 = dh(x 0 ,xi). (1.3) 

The interval (w_,w + ) C R, which can be infinite, is called the twist interval 

off)- 



Remark 1.1.2. The twist condition (1.2) states that images of verticals are 
graphs over the x-axis; see Fig. 1.1. This implies that f> can be described 
in the coordinates Xq,X\ rather than xo,yo- I n other words, for every choice 
of ^-coordinates xq and x± (corresponding to the configuration space), there 
are unique choices yo and y\ for the y-coordinates (corresponding to the 
velocities) such that the image of (xo,yo) under f> is (#i,yi). 

Remark 1.1.3. A generating function h for a twist map <f> is defined on the 
strip 

{(£, v) e R 2 1 u- < 77 - i < uj+} 




1.1 Monotone twist mappings 



3 



y 




X 



Fig. 1.1. The twist condition 



and can be extended continuously to its closure. It is unique up to additive 
constants. Equation (1.3) is equivalent to the system 



\dih(x 0l x 1 ) = -y 0 , 14 ) 

\d 2 h(x 0 ,xi) = U! 

Here, the expression di denotes the partial derivative of a function with respect 
to the i — th variable. The equivalent of the twist condition (1.2) for a generating 
function is 

did 2 h < 0. (1.5) 

Finally, a generating function satisfies the periodicity condition /i(£+l, r/+ 
1) = h(£,ri). 

Monotone twist maps are not as artificial as they might seem. They ap- 
pear in a variety of situations, often unexpected and detected only by clever 
coordinate choices. In the following, we give a few examples. The reader my 
consult 

Example 1.1.4- The simplest example is what is called an integrable twist map 
which, by definition, preserves the radial coordinate 1 . In this case, the prop- 
erty of being area-preserving implies that an integrable twist map is of the 
following form: 

<l>(xo,yo) = (xo + f(yo),yo) 

with f > 0. Then the generating function (up to additive constants) is given 

by 

1 In the context of integrable Hamiltonian systems, this means that ( x,y ) are al- 
ready the angle-action-variables. 




4 



1 Aubry-Mather theory 



h = h(x i — X’o), 

with h 1 = / _1 ; in other words, h is strictly convex. 

Example 1.1.5. In some sense the “simplest” non-integrable monotone twist 
map is the so-called standard map 

<t> : (x, y) i— > (x + y + — sin 27rx, y + — sin 27 tx) 

27T 27T 

where k > 0 is a parameter. This map has been the subject of extensive 
analytical and numerical studies. Famous pictures illustrate the transition 
from integrability (k = 0) to “chaos” (. k ~ 10). 

Example 1.1.6. A particularly interesting class of monotone twist maps comes 
from planar convex billiards; we will deal with convex billiards in Chap. 3. 
The investigation of such systems goes back to Birklroff [15] who introduced 
them as model case for nonlinear dynamical systems; for a modern survey see 
[101]- 




Fig. 1.2. The billard in a strictly convex domain 



Given a strictly convex domain 17 in the Euclidean plane with smooth 
boundary <9 17, we play the following game. Let a mass point move freely inside 
17, starting at some initial point on the boundary with some initial direction 
pointing into 17. When the “billiard ball” hits the boundary, it gets reflected 
according to the rule “angle of incidence = angle of reflection”; see Fig. 1.2. 
The billiard map associates to a pair (point on the boundary, direction), re- 
spectively ( s , i[>) = (arclength parameter divided by total length, angle with 
the tangent), the corresponding data when the points hits the boundary again. 
The lift of this map, which is then defined on R. x (0,7r), is not a monotone 
twist map. 
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However, elementary geometry shows [101] that the map preserves the 
2-form 

sin if; dtp A ds = d(— cos ip) A els. 

Hence the billiard map preserves the standard area form dx A d y in the new 
coordinates 

(. x , y ) = (s, — cos VO € M x (—1, 1). 

Moreover, if you increase the angle with the positive tangent to <912 for the 
initial direction, the point where you hit dfi again will move around dfi in 
positive direction. This means that 



dx\ 

dy 0 



> 0 , 



so the billiard map in the new coordinates does satisfy the monotone twist 
condition. 



y 




Example 1.1.7. Consider a particle moving in a periodic potential on the real 
line. According to Newton’s Second Law, the motion of the particle is deter- 
mined by the differential equation 

x(t) = V'(x(t)). 



This can be written as a Hamiltonian system 

\x{t) = d y H(x(t),y(t)) 

= -d x H(x(t),y(t)) 

with the Hamiltonian H(x,y) = y 2 / 2 — V(x). For small enough t > 0, we have 
dx(t;x(0),y(0)) d 



dy( 0) dy( 0) J Q 

/' dy(T;x(0),y(0)) 



x(t; x(0),y(0)) dr 
dr 



dy{ 0) 



> 0 . 
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Therefore the time-t-map <p * H is a monotone twist map provided t is small. 
In fact, this holds true not only for Hamiltonians of the form “kinetic energy 
+ potential energy”, but for more general Hamiltonians which are fiberwise 
convex in the second variable (corresponding to the momentum) . 

A particular case is that of a mathematical pendulum where x is the 
angle to the vertical and V'(x) = — sin 27 T£. The phase portrait in IxE, see 
Fig. 1 . 3 , shows two types of invariant curves: closed circles around the stable 
equilibrium ( “librational” circles) , and curves homotopic to the real line above 
and below the separatrices (“rotational” curves). 

Note that, by the monotone twist condition, an orbit ((a;*, j/i))iez of a 
monotone twist map 4 > is completely determined by the sequence (xj)jg z via 
the relations 

Vi = d 2 h(x i - 1 ,x i ) = -dih(xi,x i+ i). 

Similarly, an arbitrary sequence (£i)j<= z corresponds to an orbit of a monotone 
twist map </> if and only if 

c^Mfi-i, &) + dih(£i,£i+i) = 0 (1.6) 

for all i £ Z. Thus, on a formal level, orbits of a monotone twist mapping may 
be regarded as “critical points” of the discrete action “functional” 

(&)ie z 

iez 

on R z . This point of view leads to the following notion of minimal orbits. 



1.2 Minimal orbits 

Let (j) : (£0,2/0) * (£1,2/1) be a monotone twist map with generating function 

h(x o,£i). We have seen above that the </Aorbit of a point (£0,2/0) is com- 
pletely determined by the sequence (xi) of the first coordinates. Moreover, an 
arbitrary sequence (£,) corresponds to an orbit if, and only if, it satisfies the 
recursive relation ( 1 . 6 ). Loosely speaking, orbits are “critical points” of the 
action “functional” 

(£*)*e z !->• ^2 M£*>£*+i)- 
ie z 

In this section, we are interested in minima, i.e. in points which minimize the 
action. 

This, of course, makes only sense if we restrict the action of a sequence 
(£*)»ez to finite parts. In analogy to the classical Principle of Least Action, 
we define minimal orbits in such a way that they minimize the action with 
the end points held fixed. 




1.2 Minimal orbits 
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Definition 1.2.1. Let h be a generating function of a monotone twist map 
4>. A sequence (xi)jgz with G R. is called minimal if every finite segment 
minimizes the action with fixed end points, i.e., if 



i - 1 i - 1 

J2h(xi,x i+ 1) < y ‘ £i+l) 

i=k i=k 

for all finite segments £;) G R i_fe+1 with and £; = a:;. 

(1.6), eacft minimal sequence (xj)j 6 % corresponds to a (f>-orbit 
((xi,yi))i£ z; these are called minimal orbits of <f>. 

Given an orbit (xi,yi) in S 1 x ( a,b ), the twist map <j> induces a circle 
mapping on the first coordinates x*. This leads to the definition of the rotation 
number of an orbit of a monotone twist map. 

Definition 1.2.2. The rotation number of an orbit ((&», 2/i))»e z of a mono- 
tone twist map is given by 

UJ := IllTL — = lim ; 

|£| — »oo Z |i| — ^oo Z 



if this limit exists. 

Example 1.2.3. The simplest orbits for which the rotation number always ex- 
ists are periodic orbits, i.e., orbits ((&*, 2/i))iez with 

x i+q = Xi+p 

for all i G Z, where p , q are integers with q > 0. In order to have q as the 
minimal period one assumes that p and q are relatively prime. Then the 
rotation number is given by 

P 

to = -. 

q 

The questions arises whether there are orbits for a monotone twist map of 
any given rotation number in the twist interval. Actually, this is the core of 
Aubry-Mather theory, which yields an affirmative answer. The classical result 
in this context is a theorem by G.D. Birklroff [15] who proved that monotone 
twist maps possess periodic orbits for each rational rotation number in their 
twist interval. Perhaps because monotone twist maps were not that popular 
in the mid-20th century, it took 60 years to generalize Birklroff’s result to all 
rotation numbers. 

Theorem 1.2.4 (Birkhoff). Let <p be a monotone twist map with twist in- 
terval (w_,u;+), and p/q G (w_,cc+) a rational number in lowest terms. Then 
(j> possesses at least two periodic orbits with rotation number p/q. 
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Proof. The proof is a nice illustration of the use of variational methods in the 
construction of specific orbits for monotone twist maps. 

Consider the finite action functional 

m o, •■>.«&) 

2=0 

on the set of all ordered (q + l)-tuples with 

& <6 < ...<& = & + P- 

Since these tuples form a compact set, the continuous function H has a min- 
imum, corresponding to a periodic orbit of the monotone twist map <f>. What 
we need to show is that this minimum does not lie on the boundary, which 
consists of degenerate orbits of length less than q. 

Suppose that there is a periodic orbit with 

€j- 1 < = £j+i < 0+2 

for some index j: the case of more than two equal values is treated analogously. 
Then the recursive relation (1.6) yields 

£j+i) = 0 

d2h(^j,fj + i) + dih(^ j+ i,^j + 2) = 0 
Since £j = fj+i, substracting the two equations gives 

- d 2 h(t;j,£j) + dih(£j+i,£j+i) - dih(£j + i,£j +2 ) = 0 . 

This can be written as 

did 2 h(in,tj) (£j-i - $j) + d 2 d 1 h(t j+1 ,ri2) (£j+i - ^+ 2 ) = 0, 

where 771,772 are two intermediate values. But the left hand side is strictly 
negative, due to (1.6) and the assumptions, which is a contradiction. □ 

Birklioff ’s theorem is sharp in the sense that, in general, one cannot expect 
more than two periodic orbits with a given rotation number. For example, in 
the elliptical billiard, there are precisely two 2-periodic orbits, corresponding 
to the two axes of symmetry. 



1.3 The minimal action for monotone twist mappings 

Of particular importance for the dynamics of a (projection of a) monotone 
twist map <j> : S 1 x (a, b) — > S 1 x (a, b) are closed invariant curves. They fall 
into two classes: an invariant curve is either contractible or homotopically non- 
trivial. Lifted to the strip R x (a, b), this means that we consider ^-invariant 
curves which are either closed or homotopic to R. 
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Definition 1.3.1. An invariant circle of a monotone twist map <p is an em- 
bedded, homotopically nontrivial, (f-invariant curve in S 1 x ( a,b ), respectively, 
its lift to K. x (a, b) . 

Example 1.3.2. Considering the phase space 1 x K of the mathematical pen- 
dulum (see Fig. 1.3), the librational circles around the stable equilibria are 
not invariant circles according to our definition. On the other hand, the rota- 
tional curves above and below the separatrices do represent invariant circles. 
Finally, the union of all the upper, respectively lower, separatrices also form 
(non-smooth) invariant circles. 

ft turns out that invariant circles of monotone twists maps cannot take 
any form. Indeed, another classical result by G.D. Birklroff states that they 
must project injectively onto the base. More precisely, we have the following 
theorem. 

Theorem 1.3.3 (Birkhoff). Any invariant circle of a monotone twist map 
is the graph of a Lipschitz function. 

There are essentially two different proofs of this result. The original topo- 
logical approach is indicated in [15, §44] and [16, §3]; precise, and even more 
general, proofs along this line can be found in [28, 42, 51, 66, 70]. The second 
approach [94] is different and more dynamical. We give a sketch of its main 
idea here and refer to [94] for details. 

Proof ([9^]). Assume, by contradiction, that there is an invariant circle T of 
a monotone twist map <f> which is not a graph. Then we have a situation like 
that indicated in Fig. 1.4. 

Let us apply <j) once and see what happens to the area of the domain fi 0 . 
Since the preimage 1 (x^i ) is a graph in view of the monotone twist condition, 
and since (j> is area-preserving, the application of <f> pushes more area into the 
fold, i.e., the area of fl\ is bigger than that of fi 0 . 

Now iterate <fi, and consider the domains fi n for n > 1. Each application 
lets the area of fi„ grow: 

\n n \ > |c„_i| > ... > | r2i| > |fi 0 |. 

On the other hand, everything takes place in a bounded domain because T is 
an invariant curve. Therefore, we conclude that sup n |fi n | < 00 which implies 
the areas of the additional pieces tend to zero: 

lim |fi n \fi n _i|=0. 

n —> 00 

But it is easy to see that this means that r must have a point of self- 
intersection and, hence, is not embedded. 

This contradiction proves the theorem. □ 
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Let us return to the question whether there are orbits of any given rotation 
number for a monotone twist map. Theorem 1.2.4 asserts that there are always 
periodic orbits for a given rational rotation number in the twist interval. By 
taking limits of these orbits, one can construct also orbits of irrational rotation 
numbers. All of these orbits are minimal. 

Minimal orbits resemble invariant circles in the sense that they, too, project 
injectively onto the base. In other words, minimal orbits lie on Lipschitz 
graphs. Moreover, if there happens to be an invariant circle, then every orbit 
on it is minimal. 

The following theorem is the basic result in Aubry-Mather theory. The 
reader may consult [6, 34, 51, 72, 74] for more details. 

Theorem 1.3.4. A monotone twist map possesses minimal orbits for every 
rotation number in its twist interval; for rational rotation numbers there are 
always at least two periodic minimal orbits. 

Every minimal orbit lies on a Lipschitz graph over the x-axis. Moreover, 
if there exists an invariant circle then every orbit on that circle is minimal. 

Remark 1.3.5. Theorem 1.3.4 remains true if one considers the more general 
setting of a monotone twist map on an invariant annulus {(x,y) \ u_(x) < 
y < tt + (:r)} between the graphs of two functions u±\ see [72]. 

From the existence of orbits of any given rotation number, we can build a 
function which will play a central role in our discussion. Namely, consider a 
monotone twist with generating function h. Then we associate to each u> in 
the twist interval the average /i-action of some (and hence any) minimal orbit 
((xi,yi))iez having that rotation number u. 
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Definition 1.3.6. Let <fi be a monotone twist map with generating function h 
and twist interval (w_,w+). Then the minimal action of <f> is defined as the 
function a : (lv-,oj+) — > R. with 



“( y ) : = ^ h ( X i’ X i+ 1 )- 

i——N 

The minimal action can be seen as a “marked” Principle of Least Action: 
it gives the (average) action of action-minimizing orbits, together with the 
information to which topological type the corresponding minimal orbits be- 
long. We wills see in Chap. 4 how this relates to the marked length spectrum 
of a Riemannian manifold. 

Does the minimal action tell us anything about the dynamics of the under- 
lying twist map? This question is central from the dynamical systems point of 
view. It turns out that, indeed, the minmal action does contain information 
about the dynamical behaviour of the twist map. 

The following theorem lists useful analytical properties of the minimal 
action a. 

Theorem 1.3.7. Let <p be a monotone twist map, and a its minimal action. 
The the following holds true. 

1. a is strictly convex; in particular, it is continuous. 

2. a is differentiable at all irrational numbers. 

3. If u = p/q is rational, a is differentiable at p/q if and only if there is an 
^-invariant circle of rotation number p/q consisting entirely of periodic 
minimal orbits. 

4- If T u is an ((-invariant circle of rotation number to then a is differentiable 
at u) with a' (to) = f r ydx. 

Proof. Everything is well known and can be found in [72, 68], except perhaps 
for the precise value of a'(u>) in the last part. This follows from the more 
general Thm. 2.1.24 and Rem. 2.1.7 in the next section. □ 

For later purposes, we need a certain continuity property of the minimal 
action as a functional. Namely, what happens with the minimal action if we 
perturb the monotone twist map? It turns out that, at least for perturbations 
of integrable twist maps, the minimal action behaves continuously. This is 
made precise in the next proposition. 

Proposition 1.3.8. Let ho be a generating function for an integrable twist 
map such that 

h 0 (s) = c(s-'y) k + O((s- 1 ) k+1 ) 

as s — > 7 with c > 0 and k > 2. Let h be a generating function for another 
(not necessarily integrable) twist map such that 
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h(tv) = h 0 (v-0 + O((v-^-'y) k+m ) 

as r) — £ — > 7 with 2m € N \ {0} . 

TTien </ie corresponding minimal actions ao and a satisfy 

ao(w) = h 0 (u), 



as well as 

a(u) = a 0 (u) + 0((u - 7 ) fc+m ) 



as u 7. 

Proof. Let us first convince ourselves that cto = ■ This follows from the fact 

that all orbits of rotation number ui lie on the invariant circle S 1 x {(^q) _1 (w)} 
and have the same average action /i 0 (w). Hence the minimal action ao(co) is 
indeed ho(u>). 

For the continuity of the minimal action with respect to the generat- 
ing function, we will use a monotonicity argument which is standard in the 
calculus of variations; compare also [8]. Let us consider the minimal action 
a = limjv-xx) 1/2 h(xi,x i+i), where (xi) is h- minimal, i.e. , 

hfXijXi- (-i) A £i-|-l) 

for all finite sequences (£;) with the same end points. Note that the action 
of an arbitrary segment (not necessarily part of an orbit) is monotone in the 
generating function: if hi < /12 then 

*+ 1 )- 

i i 

Moreover, the minimality of a sequence (xf) is defined by a minimization 
process over all sequences (£, ), a set which does not depend on the generating 
function h. Hence, not just the action, but also the minimal action is monotone 
in the generating function. 

The monotonicity of the minimal action implies the second assertion. □ 

Later, we will apply this proposition when 7 = oj_ is the lower boundary 
point of the twist interval. Note that in this case we may have k = 3 , for 
instance, which would be forbidden if 7 were a point in the twist interval 
because then h 0 would not fulfill the generating function condition dic^o = 

-h% < 0 . 

Finally, since a is a convex function by Thm. 1 . 3 . 7 , it possesses a convex 
conjugate (or Fenchel transform ) a* defined by 

a*(I) := sup (u)I — a(w)). ( 1 . 7 ) 

U) 

Actually, a is strictly convex, so the supremum is a maximum, and a* is a 
convex, real-valued C-function with 
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(a*)'(c/(u))=u 

whenever a' (ui) exists [90, Thm. 11.13]. Flat parts of a* correspond to points 
of non-differentiability of a. 2 



See [90] for any question about smooth or non-smooth convex analysis. 



2 




2 



Mather— Mane theory 



It was well known that the theory of Aubry and Mather concerning action- 
minimizing orbits is valid only in two dimensions. For, there is a classical 
example by Hedlund [41] of a Riemannian metric on T 3 such that minimal 
geodesics exist only in three directions. Hedlund’s construction modifies the 
flat metric on T 3 in such a way that there are three directions, corresponding 
to three disjoint “highway tunnels”, along which the metric is very small, so 
that the particle can travel along these highways and gather almost no action. 
Hedlund shows that any minimal geodesic changes between the tunnels only 
finitely often. Therefore, the asymptotic directions of minimal geodesics are 
confined to the three tunnel directions. 

Hedlund’s example showed that any generalization of Aubry-Mather the- 
ory to higher dimensions could not deal with minimal orbits. Instead, Mather 
[69] developed a corresponding theory of action-minimizing invariant mea- 
sures for positive definite Lagrangian systems. Later, Mane [62] gave another 
approach using a so-called critical value. This value singled out the energy 
value at which certain dynamically relevant orbits appear. Essentially, these 
are two sides of one coin. 

In this section, we will give an introduction to the relevant notions and 
results. For further details we refer to [21, 29, 72], 



2.1 Mather’s minimal action 

The setting for Mather’s generalization of the theory of minimal orbits to 
higher dimensions are convex Lagrangian (or Hamiltonian) systems on the 
tangent (or cotangent) bundle of a compact manifold; we will restrict ourselves 
to the case of the n-dimensional torus T™ = K"/Z". 



K.F. Siburg: LNM 1844, pp. 15-35, 2004. 
© Springer- Verlag Berlin Heidelberg 2004 
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2.1.1 The minimal action for convex Lagrangians 

For the convenience of the reader, we present a quick review of the classical 
Lagrangian calculus of variations; for details we refer to [29, 33]. We denote 
by x,p the canonical coordinates on the tangent bundle TT" = T" x R". Any 
C 2 -function L : S 1 x TT" — » R, the so-called Lagrangian, gives rise to the 
Euler-Lagrange flow on TT", defined as follows. 

The action of a C' 1 -curve 7 : [a, b\ — > T" is defined as the integral 

A{i):= f L{t,j(t),^{t))dt. 

J a 

Curves that extremize the action among all curves with the same end points 
are characterized by the Euler-Lagrange equation 

= f^’ 7 ^’ 7 ^) (2- 1 ) 

for all t € [a, 6]. Equation (2.1) is equivalent to 
r ) 2 T r)T r ) 2 T 

— (tMt),i{t))l(t) = 7W,7W)7(*)- (2.2) 

If the Lagrangian satisfies the so-called Legendre condition 

, d 2 L , 

iet w* 0 ' 

then one can solve (2.2) for 7 and, therefore, define a time-dependent vector 
field X L (t,x,p) = ((x,p),(p,X£(t,x,p)) on TT" such that the solutions of 
7 (t) = X£(t, 7(t),7(t)) are precisely the curves satisfying the Euler-Lagrange 
equation (2.1). The vector field X l is called the Euler-Lagrange vector field, 
and its flow is called the Euler-Lagrange flow p>L- It turns out that pl is C 1 , 
even if L is only C 2 . 

Definition 2.1.1. A convex Lagrangian is a C 2 -function 

L : S 1 x TT" -► R 

such that the following conditions hold. 

1. Restricted to every fiber {t} x T x T", L is strictly convex; this means that 
L has fiberwise positive definite Hessian: 
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2. L has fiberwise superlinear growth (with respect to some, and hence any, 
Riemannian metric on T n ); this means that 



lim 

|p|— >oo 



L{x,p) 

\p\ 



= 00 



uniformly in x. 

3. The Euler-Lagrange flow ipi, is complete, i.e., its solutions exist for all 
times. 



Example 2.1.2. A prime example of a flow generated by a convex Lagrangian 
is the geodesic flow on T" with respect to some Riemannian metric, where one 
considers the free motion of a particle on T”. The Lagrangian is then given 

by 

L{x,p) = i \p\ 2 x . 

If one adds a potential V on T", the Lagrangian changes to 

L(x,p) = i \p\ 2 x - V(x). 

Remark 2.1.3. A Lagrangian is by no means uniquely defined by the Euler- 
Lagrange flow. Indeed, if L generates the flow ipL, then also the new La- 
grangian 

L(x,p ) - v x [p), 

where v is any closed 1-form on T", generates the same (fiL- 

This can be seen as follows. The actions of a curve 7 with respect to L 
and L — v differ by the term v. Since v is closed, Stokes’ Theorem im- 
plies that this term does not depend on the curve 7 (in the same homotopy 
class). Therefore the actions differ only by some additive constant, and so the 
extremal curves are the same. 

Note that for convex L , the new Lagrangian L v is also convex. 



Let L be a convex Lagrangian. In the following, we will not deal with 
orbits of the Euler-Lagrange flow ipL, but rather with invariant probability 
measures. To do so, we denote by Ml the set of ipL -invariant probability 
measures on TT™. For p € Ml we call 




£ M U {+ 00 } 



its action. To each p G Ml, one associates the linear functional 



H\ T”,K) -> R. 



I — > 



/ 



v dp 



where we view a 1 -form v as a function on T T" that is linear on the fibers. 
By duality, there is a unique class p(p) G Hi (T" , R) such that 
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J vdp = ([v],p(fj,)) (2.3) 

for all [v] £ H 1 ( T",R). 

Definition 2.1.4. Let p £ Ml be an invariant measure for a convex La- 
grangian L. Then the class p(p) € .ffi(T",R), defined by (2.3), is called the 
rotation vector of p. 

Remark 2.1.5. The rotation vector of an invariant measure is related to 
Schwartzman’s asymptotic cycles [91]; see [21]. 

In analogy to Aubry-Mather theory in two dimensions, we want to min- 
imize the action of all invariant measures having the same rotation vector. 
Although the tangent bundle of T" is not compact, this can be dealt with 
by taking its one point compactification, adding a point at infinity; see [69]. 
Then Ml becomes compact with respect to the vague (weak*) topology [14], 
and we actually can minimize the action over the set of invariant probability 
measures having a given rotation vector. 

Definition 2.1.6. Let L be a convex Lagrangian. Then the function 

a : H 1 (T n , R) -> R 

h i— > min{A(^) | p € Ml,p(p) = h} 

is called the minimal action of L. 

Any invariant measure p £ Ml realizing this minimum, i. e. with A(p) = 
a(p(p)), is called a minimal measure. For a fixed rotation vector h £ Hi (T n , R), 
the set of all minimal measures with p(p) = h is denoted by Mh- 

Remark 2.1.7. In the case of one degree of freedom (n = 1), the theory of 
Mather Mane reproduces the discrete Aubry-Mather theory from Chap. 1. 
To see this, one uses the result by Moser [78] that every monotone twist 
map on the cylinder is the time-l-map of a convex Lagrangian; see also [93]. 
Then it is shown in [67] that the minimal action a(p(p)) in the continuous 
setting considered here is, perhaps after adding a constant, the same as the 
minimal action a (to) in the discrete framework of Aubry-Mather theory where 
p(p') = co. Hence we need not distinguish between the two. 

Remark 2.1.8. The relation between minimizing measures and globally mini- 
mizing orbits is quite delicate, and we refer to [21] for details. We mentioned 
Hedlund’s example [41] showing that minimal orbits for an arbitrary rotation 
vector need not always exist. At least, every trajectory that lies in the union 
of all supports of minimal measures in Mh minimizes the action among all 
curves in the universal cover R" with the same end points [69, Prop. 3]. The 
dynamics on the set of minimizing trajectories is not limited to any particular 
behaviour — it can be as complicated as that of any vector field on the base 
manifold [60]. 
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Let us consider the minimal action a. Recall from Thm. 1.3.7 that, in the 
two-dimensional discrete setting of Aubry-Mather theory, the minimal action 
is a strictly convex function. We want to prove a similar result for the higher 
dimensional case. 

Proposition 2.1.9. The minimal action a : Hi( T ra ,K) — > R. is a convex, 
superlinear function. 

Proof. Let ft. 1,/12 £ J?i(T”, R) and A £ [0,1]. Choose minimal measures 
Pi,P 2 £ Ml such that p(pi) = hi. Then the convex combination 

p := A/zi + (1 — A)/i2 

lies in M L and has rotation vector p(p) = Xhi + (1 — A)/i 2 - Since both pi and 
P 2 are minimal, we conclude that 

a(Xh\ + (1 — A) 2 ) ^ A(/r) = Xa(hi) + (1 — A)o;(/i2), 

which proves the convexity of a. 

As for the superlinearity, we refer to [69] or [29] [Thm. 4.4.5]. □ 

Remark 2.1.10. In contrast to the two-dimensional case, the function a need 
not be strictly convex. 

As a convex function, a possesses a convex conjugate 

a* : H 1 { r\M)->]R (2.4) 

c 1 — > sup ((c, h) — 01 (h)) (2.5) 

heH 1 

Since a is superlinear, the supremum is a maximum and attained at h c £ 
Hi(T n ,M) if, and only if, 



a(h) > a(h c ) + (c,h — h c ) 

for all h , in other words, if c is a subgradient of a at h c ; compare, for instance, 
[90, 29]. We arrive at the following equivalent formulations for the minimality 
of a measure p: 

• there exists a homology class h £ H i(T",K), namely the rotation vector 
p(p), such that p minimizes the action J L dp amongst all measures in 
A 1l with rotation vector h\ 

• there exists a cohomology class c = [v\ £ H 1 ( T",R), namely any subgra- 
dient of a at p(p), such that p minimizes f L — v dp amongst all measures 
in Ml- 

Note that L — v is again a convex Lagrangrian and generates the same flow 
as L because v is closed. Therefore, Ml-v = Ml', see Rem. 2.1.3. 

Let us continue with the idea to prove results, analogous to those in Aubry- 
Mather theory, in the more general setting of Mather’s theory of minimal 
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measures. Recall that Thm. 1.3.4 stated that minimal orbits of monotone 
twist maps always lie on Lipschitz graphs. Thus, one is lead to the conjecture 
that the supports of minimal measures (corresponding to minimal orbits) 
should lie on Lipschitz graphs over T" (seen as the zero section in TT"). 

In fact, this conjecture is true. The following is Mather’s so-called Lipschitz 
Graph Theorem from [69]; see also [21]. 

Theorem 2.1.11. For every h £ Hi( T",R), the union of the supports of all 
minimal measures in A ih lies on a Lipschitz graph over T". Moreover, the 
Lipschitz constant depends only on the Lagrangian L and not on the rotation 
vector h. 

Important dynamical objects for twist maps are invariant circles; in higher 
dimensions, the corresponding objects are invariant tori. We know from 
Thm. 1.3.4 that orbits on invariant circles are automatically minimal. What 
is the corresponding result in higher dimensions? We point out that invariant 
tori of convex Lagrangian systems are only shown to be graphs under cer- 
tain assumptions on their dynamics; see [10] for a generalization of Birklroff ’s 
Theorem 1.3.3 to higher dimensions. 

In order to deal with invariant tori, it is convenient to reformulate every- 
thing in the Hamiltonian, rather than in the Lagrangian, framework. Given a 
convex Lagrangian L : S 1 x T T™ — > R, the so-called Legendre transformation 

t : S 1 x TT n -> S 1 x T* TP 

( 2 . 6 ) 

(t,x,p) (t, x, y := d p L ) 

is a diffeomorphism between the tangent and the cotangent bundle. It yields 
the convex Hamiltonian H : S 1 x T* T" — > R defined by 

H(t,x,y ) := (y,p) - L(t,x,p) \ P =( dpL) -i( y )- 

The Hamiltonian H gives rise to the Hamiltonian flow pn on the cotangent 
bundle via the Hamiltonian equations, written in local coordinates as 



x(t) = d y H(t,x(t),y(t )) 
y(t) = -d x H(t,x(t),y(t)) 



(2.7) 



Then the Legendre transformation provides a conjugation between the Hamil- 
tonian flow (fix on T* T" and the Euler-Lagrange flow ipi, on T T". We refer 
to [33, 21] for more details. 

Given a Hamiltonian flow pu on T* T n , we denote its time-t-map by 



<Ph 



rp*rjpr) 






This yields a one-to-one correspondence between ^-invariant probability 
measures and ph~ or (^-invariant ones. For simplicity, we do not introduce 
three different notations but write p for any of those. Likewise, we define the 
minimal action associated to a convex Hamiltonian H to be that associated to 
L and write a in either case. We say that a Ph~ or ^5^-invariant probability 
measure is minimal if its Pl -invariant counterpart is. 
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2.1.2 A bit of symplectic geometry 

The Hamiltonian viewpoint is the viewpoint of symplectic geometry. Let us 
recall a few notions; see [73] for a comprehensive introduction to symplectic 
geometry. 

Definition 2.1.12. A symplectic form u> on a manifold M is a closed non- 
degenerate 2-form. A symplectic manifold is a manifold M, equipped 

with a symplectic form u>. 

Example 2.1.13. The 2n-dimensional Euclidean space R 2 ", together with the 
so-called canonical symplectic form 



o>o : = dy A d.T = dpi A da c*, 

i—1 

is called the standard symplectic space. Note that the dimension of a symplectic 
manifold must always be even in view of the nondegeneracy condition on the 
2-form uj. 

Example 2.1.14. An important example of a symplectic manifold is the cotan- 
gent bundle T*X of an n-dimensional manifold X. It carries a canonical 
symplectic form w = dA that is not just closed but even exact. Here, the 
1-form A is the so-called Liouville form which, in local coordinates, is given 
by 

n 

A := y da; = ^yjda;;. 

i=i 

This local definition admits a global interpretation as follows. Let 

9 : T*X -> X 

be the canonical projection, and £ G T^ x y ^T*X. Then 

A(x,v)(0 = (**ifc)(0- 

Of particular interest in symplectic geometry are submanifolds A C M of 
s symplectic manifold ( M,u > ) on which the symplectic form vanishes: 

w\ta = 0 . 

Such submanifolds are called isotropic. It follows from the nondegeneracy of 
uj that dim A < 1/2 dim M for isotropic submanifolds A. 

Definition 2.1.15. A Lagrangian submanifold A of a symplectic manifold 
( M , uj) is an isotropic manifold of maximal dimension; in other words, we 
have 

dim A = — dim M and uj\ta = 0. 
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Example 2.1.16. In the standard symplectic space, the submanifold {(x,y) € 
R 2n | y = 0} is a Lagrangian submanifold, whereas {(x,y) G K 2 " | x = 0} is 
not. 

Example 2.1.17. Let v be a 1-form on some manifold X. Then the graph 

gri/ := {(x,u x ) | x G X} 

is a Lagrangian submanifold of (T*X, dA) if, and only if, the 1-form v is 
closed. Such a Lagrangian manifold, which projects injectively onto the base, 
is called a Lagrangian graph or Lagrangian section. 

In our case where M = T*T™, any Lagrangian submanifold that is diffeo- 
morphic to T™ is called a Lagrangian torus. For instance, if n = 1, any circle 
on the cylinder is a Lagrangian torus (or circle, rather) . 

We want to define Hamiltonian flows on symplectic manifolds. To do so, 
let H : S 1 x M — > R. be a time-periodic Hamiltonian on some symplectic 
manifold ( M,u> ), and denote by H t : M — > R. the function for fixed t. 

Definition 2.1.18. The Hamiltonian vector field Xh on M associated to a 
Hamiltonian H is defined by 



ix H u = -d H t , 

where ix H u '■= u(Xh, •) is the usual contraction of a form by a vector field. 

Example 2.1.19. If = (R 2 ",wo) is the standard symplectic space then 

the Hamiltonian vector field is given by 

X H (x,y) = JVH t (x,y), 

where J is the 2 n x 2?r-matrix 




In other words, we arrive at our familiar system (2.7). 

The invariance group of a symplectic manifold consists of all diffeomor- 
phisms that leave the symplectic form invariant. 

Definition 2.1.20. A map <f> : M — > M of a symplectic manifold is 

called symplectic if it preserves the symplectic form u>: 

(j)*LO = UJ. 

Example 2.1.21. Certainly, the identity is symplectic. More generally, every 
time-f-map tp ^ of a Hamiltonian flow is symplectic since, by Cartan’s formula 
for the Lie derivative, we have 

^(<^)*w = L Xh u) = d ix H u> + ix H duj = d(-d H t ) = 0. 
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Example 2.1.22. On a cotangent bundle T*X with coordinates (x,y) and 
canonical symplectic form (see Ex. 2.1.14), we have the symplectic shift map- 
ping 

{x,y) (x,y- v) 

where v is some closed 1-form on X. 



2.1.3 Invariant tori and the minimal action 

Let us return to our original setting. We know from Thm. 1.3.4 that invariant 
circles of monotone twist maps carry minimal orbits. In higher dimensions, 
a similar statement is true. Namely, let (j> = ipjj be generated by a convex 
Hamiltonian on S 1 x T*T™, and suppose that f> possesses an invariant La- 
grangian torus A which is a graph. This situation occurs, for instance, in 
KAM-theory where one considers small perturbations of convex, completely 
integrable Hamiltonian systems. 

Definition 2.1.23. Consider a cotangent bundle 6 : T*X — > X with its 
canonical symplectic form to = dA. We denote by C the class of all Lagrangian 
submanifolds ofT*X which are Lagrangian isotopic to the zero section O. 

Given A £ C, the natural projection 9\a '■ A — ■> X induces an isomorphism 
between the cohomology groups H 1 (X, R) and H 1 (A,K). The preimage a a € 
H x {X, HR) of [A | a] € i/ 1 (A, R) under this isomorphism is called the Liouville 
class of A. 

The next theorem, firstly, says that A consists of supports of minimal 
measures and, secondly, shows that the Liouville class of A is a subgradient of 
the minimal action. Recall that a vector v £ R" is a subgradient of a function 
/ : R n — > R at cc G R n if 



f(y) > f(x) + {v, y - x) 

for all y £ R”. If we have a strict inequality for all we say that v is a 

subgradient with only one point of tangency. For instance, if / is differentiable 
at x then, of course, its gradient V/(x) is its unique subgradient at x. See 
[90] for more details. 

Theorem 2.1.24. Let 4> = ip\j be generated by a convex Hamiltonian H on 
S 1 x T* T". Suppose that <f) possesses an invariant Lagrangian torus A in 
(T*T™, d\) such that A is homologous to the zero section and <j >\ a is conjugated 
to a translation on T" by some fixed vector p. 

Then every (f-invariant. probability measure with support, in A is minimal, 
and a a € U 1 (T",R) is a subgradient of the minimal action a of H at p with 
only one point of tangency. Vice versa, every minimal measure of rotation 
vector p has support, in A. 
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We point out that an observation by Herman [43, Prop. 3.2] shows that 
the condition on A being Lagrangian can be dropped if the vector (1 , p) is 
rationally independent, e.g. for invariant KAM-tori; in this case the minimal 
measure supported on A is unique. 

Proof. We proceed in three steps and reduce each to the previous one. First of 
all, by a higher-dimensional version of BirklrofT’s Theorem [10] the tori ip 1 ^ (A) 
are graphs for all t £ [0, 1]; for our assumption on p\a says, in particular, that 
p \ a preserves a measure which is positive on open sets (cf. [10, Prop. 1.2. (ii)]). 
Note that, as a Lagrangian graph, p* H (A) is the graph of a closed 1-form up. 
by invariance, u 0 = u\. 

Case 1 : Our starting point is the simplest possible, where A = O is the 
zero section and remains invariant under the flow, i.e. 

v t = 0 

for all t. Then 

0 = y(t) = d p L(t,x(t),x(t)) 

and 

0 = -^<9 p L(t,x(t),x(t)) = d x L(t,x(t),x(t)) 

for all orbits starting (and hence lying) on O. Note that £~ 1 0, the preimage 
of O under the Legendre transformation, will depend on t unless x(t) = p for 
all t. 

In any case, we have VLt \t-^o = 0 which, by convexity of L, implies that 
L t (x,p ) = min L t 4=> (x,p) £ £~ 1 0. 

Consequently, an invariant measure p £ Ml is globally minimizing if, and 
only if, its support lies in £~ 1 0. Since all orbits in O have rotation vector p 
we see that for any such p we have 

J L dp = J min L t dt = a(p). 

In addition, a{h) > a(p) if h ^ p, so 0 = ao is a subgradient of a at p with 
only one point of tangency. 

Case 2: Next we consider the case when A is still the zero section but does 
not stay invariant under the flow; more precisely, we assume that 



v t = dS t 



for some function St on T" with v o = iq. (St is a generating function in the 
simplest case where the Lagrangian is a graph.) 

We define the new Hamiltonian 



K(t, x, y) = d t S t (x ) + H(t,x,y + dS t {x)). 
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It is convex, and if we write <^(a :,y) = (x(t),y(t)) the transformation law of 
Hamiltonian vector fields yields 

<pk( x > y) = <Pd t st ° vh( x ^ y) = (*(*)> y(*) ~ ds t (x(t))). 

Therefore, p X K = p x H = p and p t K (x,0) = (x(t),0). 

Now we are in the first case with H replaced by K. This changes the 
(minimal) action only by an additive constant [67], and the same conclusions 
hold as before. 

Case 3: In the general case, we apply the symplectic shift 

(x,y) ^ -> (x,y- n 0 ). 

This maps A onto the zero section O and the 1-form A onto A = A — A|yi; the 
new flow p*~ maps O onto the graph of = vt ~ Vo- A generating function 
for is given by 

St(x)= [ -H{s,p s ~ o(p t ~)- 1 (x,u t {x))) ds. 

Jo 

Thus we are back in Case 2, but this time with with a different 1-form A 
instead of A. It is shown in [67] that the actions behave like 

My) = My) + - A], p(M) = My) - {a a, p{y)). 

From Case 2 we know that every /i with support in A minimizes A(/u) among 
all measures with p{y) = p- But under this constraint the correction term 
(a a, p) is a mere constant, so p, minimizes A(p), too. Moreover, 0 is a sub- 
gradient of a(h) = a(h) — ( aA , h) at p with only one point of tangency. That 
means that 

a(h) > a(p) + (aA, h — p) 

with equality only for h = p, so aA is a subgradient of a at p with only one 
point of tangency. 

This finishes the proof of the theorem. □ 

Corollary 2.1.25. If, under the assumptions of Theorem 2.1.24, the invari- 
ant torus A is invariant under the flow of H (and not just its time 1-map) 
then 

a* (a a) = [ \a dt. 

Jo 

Proof. Theorem 2.1.24 implies that a*(au) = (aA, p) — Mp) = / — L dp 

for every invariant measure p supported in A. Since A is invariant under the 
Hamiltonian flow, the function H(t, •) is constant on A with H = A|/i — L. □ 

Remark 2.1.26. If H is autonomous, Corollary 2.1.25 determines the energy 
level of invariant tori with given a*(ayi); compare [25]. 
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2.2 Mane’s critical value 



Another approach to a generalization of Aubry-Mather theory to higher di- 
mensions was suggested by Mane [62] . Its main idea is to single out a certain 
energy level at which a significant change of the dynamical behaviour takes 
place. This produces a “critical” energy value for each convex Lagrangian. 
It turns out that this value is the minimum of the actions of all invariant 
measures in Ml, and that one can recover Mather’s minimal action from it 
(and vice versa). 



2.2.1 The critical value for convex Lagrangians 

Let L : XT" -> R be a time-independent convex Lagrangian on the tangent 
bundle of the n-torus. Let t : T T" — » T*T" be the Legendre transformation, 
and H : T* T" — >■ R the Hamiltonian corresponding to L. The push-forward 
of the Euler-Lagrange flow ifL on XT" by the Legendre transformation is 
the Hamiltonian flow ipn on X*T" with respect to the canonical symplectic 
structure on X*T"; see Sect. 2.1.2. The energy of L is the function E : XT" — > 
R defined by 



E(x,p) := (—{x,p),p} - L(x,p) = H(£(x,p)). 

It is a first integral of the Euler-Lagrange flow ifL- 

Recall that a curve 7 : [a, b } — > T" is called absolutely continuous if for 
every e > 0 there exists 6 > 0 so that for each finite collection of pair- 
wise disjoint open intervals in [a,b] of total length less than <5 one has 

JAdist(7(tj),7(sj)) < e; here dist is any Riemannian distance on T". As be- 
fore, the action of an absolutely continuous curve 7 : [a, b] — > T™ is defined 
by 

a l{ 7) := f L(^(t),i{t))dt. 

J a 

We keep the subscript L in order to distinguish between the actions for dif- 
ferent Lagrangians. 

Given two points X\, X2 G T" and some T > 0, denote by Ct{ 27, £2) the set 
of absolutely continuous curves 7 : [0, T } — » T™ with 7(0) = X\ and 7 (T) = x 2. 
For each fceR, we define 

&k(xi , x 2 ;T) := inf{A L+fc (7) | 7 e C T (xi,x 2 )} 

as the infimum of the ( L + fc)-actions over all curves connecting 37 and x 2 in 
time T. 

Definition 2.2.1. Let L : 7T" — > R be a convex Lagrangian. Then the action 
potential ^ : T" x T" — > R U {—00} of L is defined as 
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<P k : T" x T” — > M U {- 00 } 

(xi,x 2 ) >-> mf 4> k (xi,x 2 -,T) 

The critical value of L is given by 

c(L) := inf{fc € M | <& k (x,x) > —00 for some x € T"}. 



x 




Remark 2.2.2. Observe that if there is some closed curve 7 with AL+ k {l) < 0 
then immediately l P k (x,y) = —00 for all x,y. This follows by considering the 
curve going from x to some point on 7 , going around as many times as one 
wishes (gathering as much negative action as one wants), and finally going to 
y\ see Fig. 2.1. Thus, we could replace the word “some” in the definition of 
c(L) by “all”. Since the action potential is monotone in k, we then have 

c(L) = sup{fc € M | there is a closed curve 7 with A L+k ( 7 ) < 0}, (2.8) 

which gives another description of the critical value. 

Remark 2.2.3. We will explain why c(L) is a real number. Think of some 
Lagrangian, and pick a point x GT n . Consider the infimum of the actions of 
all closed curves through x. Since the time interval is free, you will get —00 for 
the infimum as soon as you have just one closed curve with negative action. 

Now let L be a fixed Lagrangian. We want to see what happens if we 
shift L by some constant k. If k < — min x L{x , 0) then L(x, 0) + k is negative 
at some point x, we can choose the constant curve at x , and end up with 
<Pk(x,x) = — 00 . On the other hand, the fact that L is convex implies that L 
is bounded from below. Therefore, if k > — min L then L + k is positive, and 
we must have T>k{xi,X 2 \ T) > 0 for all 27 , x 2 , T. This shows that 

c{L) < 00 



is a real number. 
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The critical value can in fact be characterized in a variety of ways [60, 20, 
22, 23]. Each of these characterizations gives new insight into the geometry 
or the dynamics of the given Lagrangian system. In the following, we will 
explain the relation between the critical value and the minimal action defined 
in Sect. 2.1. 

Let Ml be the set of invariant probability measures on T T". The next 
result, due to Mane [60], states that the critical value of a convex Lagrangian is 
equal to the minmal action of all measures in Ml, regardless of their rotation 
vector. 

Proposition 2.2.4. 



c(L) = - xmn{A L (fL) | /i G Ml}- 

Proof. First of all, one can show that 

min{zli(/x) | p G Ml} = min{2lz / (^ 7 ) | 7 abs. cont. curve} 

where p 7 is the measure equally distributed along some absolutely continuous 
curve 7; see [21]. So we will prove that 

—c(L) = min{Ai(^ 7 ) | 7 abs. cont. curve}. 

For any curve 7, we have d t+c ( t )(/i 7 ) > 0 by definition of c(L). Therefore, 

— c(L) < min{^4 i (^ 7 ) | 7 abs. cont. curve}. 

To prove the reversed inequality, we observe that, whenever k < c(L), there 
exists a curve 7 with AL+k(fij) < 0, which implies 

— k > min{Ai(/Lt 7 ) | 7 abs. cont. curve}. 

Now let k tend to c{L). □ 

Remark 2.2.5. The fact that there is an invariant measure /i with Al(p) = 
mi{A L (p) \ fie Ml} follows from the compactness of Ml as in Sect. 2.1.1. 

Recall that 

a*{c) = max ((c, h) — a(h)) 

hdzHi 

is the convex conjugate of the minimal action a : Hi(T",M) — > R defined in 
Def. 2.1.6. Therefore, 

a*(0) = -minmin {A L (li) \ p{fi) = h} = - mm{A L (fi) \ fi G Ml}, 



which yields the following description of the (convex conjugate of the) minimal 
action in terms of the critical value. 
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Corollary 2.2.6. For every closed 1-form v on T", we have 

c{L — v) = a*(\v\). 

Thus, at least in the case of the torus, the theories of the minimal action 
and the critical value are equivalent. 

It turns out the critical value c(X) of a convex Lagrangian can be, in fact, 
recovered also from the Hamiltonian H , as the following result [22] shows. 
Namely, we have 



c(L) = inf max H (x , du(x)) . (2-9) 

■ueC~(T’\R)a;eT» 

In other words, the critical value is a minimax value of H over all exact 
Lagrangian graphs. We will give a purely symplectic description of the critical 
value in Sect. 6.3.2. 

Let L : T T" — > R. be a convex Lagrangian, and M l denote the set of 
invariant probability measures on T T". 

Definition 2.2.7. A measure /io € Ml is called globally minimizing if it 
minimizes the action amongst all invariant measures, i.e., if 

A L {p o) = min {A L {p) I M € Ml}- 

The Mather set in XT™ is defined as the closure of the union of the supports 
of globally minimizing measures: 

M := U {supp(p) | /i globally minimizing}. 

Note that a globally minimizing measure must have zero rotation vector. 
Therefore, in view of Thm. 2.1.11, the set M is a Lipschitz graph with respect 
to the canonical projection 



t : TT n -> TP. 



We call the set 

M := t(M) 

the projected Mather set. It is known [25] that M is contained in the energy 
level X -1 (c(L)). Finally, we define the Mather set in T* T" as the image of 
M under the Legendre transform: 

M * 



2.2.2 Weak KAM solutions 

In the study of the dynamics of a Lagrangian system on XT", a particular 
role is played by invariant tori. Since we consider time-independent convex 
Lagrangians, we know that the energy E[x,p) on XT" (or the Hamiltonian 
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H(x,y) on T*T") is an integral of the corresponding flow. If an invariant 
torus is a graph of some closed 1-form u, then v must satisfy the equation 
H(x, v x ) = k = const. Actually, we can restrict ourselves to exact 1-form by 
considering L — v instead of L ; compare Rem. 2.1.3. Finding a smooth exact 
invariant torus gr du, where / : T" — > R. is some smooth function, is then 
equivalent to finding a smooth solution of the (time-independent) Hamilton- 
Jacobi equation 

H(x,du(x)) = k (2-10) 

Now, smooth solutions of the Hamilton- Jacobi equation will, in general, 
not always exist. There is a general theory, developed by Fatlri [29], that deals 
with Lipschitz continuous solutions of (2.10) and gives another approach to 
the critical value. Recall from (2.9) that the critical value is characterized as 
the infimum of energy values k such that the sublevel set {H < k } contains a 
smooth solution of (2.10). Moreover, the only energy level that might support 
a C 1+Llp -solution of (2.10) is the level H~ l (c(L)) where c(L) is the critical 
value. 

Example 2.2.8. It might be instructive to consider the simplest nontrivial case. 
Let L : T S 1 -^Rbe the Lagrangian of the mathematical pendulum: 

L{x,p) = i|p| 2 -cos(27ra:). 

The corresponding Hamiltonian is given by 

H{x, y) = i|i/| 2 - cos(27re) 

whose level sets are shown in Fig. 2.2. Since we are in one degree of freedom, 
the level sets of H consist of solutions of the Hamiltonian flow, and we end 
up with the phase portrait of the pendulum (Fig. 1.3). 




Note that the energy of the separatrices is 1. For energies below 1, each 
level set is a closed C 00 - curve around the point (0, 0); thus, no level set con- 
tains even a graph. For energies above 1, each level set consists of two C°°- 
curves that go around the cylinder, one above and one below the separatrices; 
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but none of these curves is exact, i.e., the graph of an exact 1-form. The 
energy level H~ 1 ( 1), however, is the union of two graphs of differentials of 
C'l+L'P-finictions u± on S 1 , namely the upper and lower separatrix. 

Therefore, the critical value is 1. In fact, c(L) = max V for any Lagrangian 
L(x,p) = 1/2 \p\1 — V(x) on TT™. 

Given a continuous function u : T" — > K, we write 

u -< L + c 

whenever u(x) — u(y) < <& c (y,x) for all x,y € T™. Here, <P C is the action 
potential for the critical value 



c := c(L). 

Remark 2.2.9. Fathi showed that a function u satisfies u -< L + c if, and only 
if, it is Lipschitz continuous and fulfills the inequality 

H(x, du(x)) < c 

for almost every x € T"; see [29]. In other words, for Lipschitz continuous 
functions u, the condition u -< L + c is equivalent to u being a subsolution of 
the Hamilton- Jacobi equation (2.10). 

Note that, by Rademacher’s theorem, Lipschitz functions are differentiable 
almost everywhere. 

Definition 2.2.10. We say that a continuous function u+ : T" — » R is a 
positive weak KAM solution if u+ satisfies the following two conditions: 

1. u+ -< L + c 

2. for all x € T" there exists an absolutely continuous curve 7^ : [0, 00) — > T™ 
such that 7® (0) = x and 

u+(7+(*)) ~u+{x) = [ (L + c)(7® (s), 7“ (s)) ds 
Jo 

for all t > 0. 

Similarly, we say that a continuous function U- : T" — > M is a negative 
weak KAM solution if u- satisfies the following two conditions: 

1. U- -< L + c 

2. for all x € T™, there exists an absolutely continuous curve 7® : (— 00, 0] — » 
T" such that 7® (0) = x and 



u-(x) - u_(7*(-t)) 



(i + c)(7K«),7^W)ds 



for all t > 0. 
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We denote by S± the set of all positive (respectively, negative) weak KAM 
solutions. A pair of functions (u-,u+) is said to be conjugate if u± G S± and 
U- = u+ on the projected Mather set A4. 

Fathi’s Weak KAM-Theorem asserts that positive and negative weak 
KAM solutions always exist [29]. Moreover, at any point £o of differentia- 
bility of a weak KAM solution u, Conditions 1 and 2 of the above definition 
imply that u satisfies 

H(x 0 ,du(x 0 )) = c. 

In fact, the points xo of differentiablity of u+ (respectively, u_) are precisely 
those for which the curve 7“° (resp. 7* 0 ) is unique. The following result [29, 
Thm. 5.1.2] states that any function u with u -< L + c can be squeezed between 
a (unique) pair of conjugate functions. 

Theorem 2.2.11. If u : T™ — > R. is a function such that u -< L + c, then there 
exists a unique pair of conjugate functions (u-,u+) such that 

u+ <u< U-. 

Finally, we will need the notion of Peierls barrier which goes back to 
Mather [71]. 

Definition 2.2.12. The Peierls barrier h : T” x T" — > R. is defined by 
h(x, y) := liminf 4> c (x,y;T). 

T — >-oo 

The function h is Lipschitz continuous and, by definition, satisfies 

h{x, y) > $ c (x,y). 

Moreover, the Peierls barrier can be written as 

h(x,y) = max (u-(y) — u+(x)) (2.11) 

(u_,u + ) 

where the maximum is taken over all pairs (u-,u+) of conjugate functions; 
see [29]. 

2.2.3 The Aubry set 

By definition, two conjugate functions u± coincide on the projected Mather 
set M. . In general, however, there is a bigger set with this property. To define 
this set we set 



2: («_,«+) := {x G T" | u-(x) = u+(x)}. 
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Definition 2.2.13. The projected Aubry set A is the set of points in T" at 
which all pairs of conjugate functions coincide: 

A . |^| 

(u_,u + ) 

where the intersection is taken over all pairs of conjugate functions. 

It follows from the definition of conjugacy that A i C A. 

In order to define the Aubry set in T* T", we note that the functions u_ 
and are differentiable at every point x € T ( u _ with the same derivative. 
Moreover, the map 

l ( u_,u + ) du-(x) = du+(x ) € T* T" 

is Lipsclritz continuous [29, Thm. 5.2.2]. That map defines a set 

l(u-,u + ) C T*T n 

that projects injectively onto I( u _ tU+ ) and contains the Mather set. The Aubry 
set in T* T” is defined as 

A* := f) i iu _,u + ) 



where, again, the intersection is taken over all pairs (u—,u+) of conjugate 
functions. As the notation suggests, one can prove that 

A = 0(A*) 

where 9 : T*T n — > T" is the canonical projection. As usual, we denote the 
preimage of A* under the Legendre transform by A and call it the Aubry set 
in IT”. 

The sets M. and A are compact and invariant under the Euler-Lagrange 
flow (j>t- It turns out that the Aubry set consists of a distinguished kind of 
orbits. To make this precise, we give the following definition due to Mane. 
Recall that c = c(L) denotes the critical value of the Lagrangian L. 

Definition 2.2.14. We say that an absolutely continuous curve 7 : [a, b] — > 
T” is semistatic if 

Al+ c ( t ) = &c(7(a),j(b)). 

An absolutely continuous curve on an infinite interval is called semistatic if 
it is semistatic on every finite interval. 

Semistatic curves are solutions of the Euler -Lagrange equation because 
of their minimizing properties. It is not hard to check that semistatic curves 
have energy precisely c. 
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Definition 2.2.15. An absolutely continuous curve 7 : [a, b] — > T" is called 
static if it is semistatic and satisfies 

7(b)) + $ c {l(b), 7(a)) = 0. 

An absolutely continuous curve on an infinite interval is called static if it is 
static on every finite interval. 

In fact, one does not even need to require that a static curve be semistatic. 
If one defines static curves by the condition that Al + c { 7) = —^(7(6), 7(a)) 
then it follows that a static curve is semistatic [21]. 

The following proposition gives a useful characterization of the Aubry set. 
Its proof is well known to experts; nevertheless, we include it for the sake of 
completeness. 

Proposition 2.2.16. The Aubry set A consists precisely of those orbits whose 
projections to T" are static curves. 

Proof. Take (x, v ) G A. We want to show that the curve 

l(t) ■= T(</> t (x,v)) 

is a static curve. By the definition of A and [29, Thm. 5.2.2] we have for any 
pair (u-,u+) of conjugate functions that 

u+il (t)) - u-(y(s)) = u + (y(t)) - u+(y(s)) = A L+c ( 7 | [S)t] ) 

for all s < t. Using (2.11) we can choose a pair (u-,u+) of conjugate functions 
for which the Peierls barrier h satisfies 

MtW,7(s)) = u-(y (s)) - u+(y(t)). 

Therefore, we can estimate 

A L+c{y\[ s ,t])+^c(y{t),y{s)) < u + (y(t))-u-(y(s))+h(y(t),y(s)) = 0. (2.12) 
It is easy to show that <P C satisfies the triangle inequality 
&c(x, V ) < $c(x, z ) + $ c (z, y) 
as well as <d> c {x,x) = 0 for all x G T™. Hence we have 

0 = &c(7(s),y(s)) 

< &c(y(s),y(t)) + & c {y{t),y{s)) 

< A L+c ( 7 | [Sjt] ) + <f> c ( 7 (t), 7 (s)) 

< 0 



in view of (2.12). This implies that 7 is a static curve. 
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Suppose now that 7 : R. — > T" is a static curve. Then 7 is a semistatic curve 
with energy c, and given s < t and e > 0, there exists a curve 7 connecting 
7 (t) to 7(s) such that 



Al+c (7I [s,t]) + A L+C {y) < e - 

Looking at the loop formed by 7|[ s ,t] and 7, we conclude that h(y(t),'y(t)) < 0. 
But h( x,x) > <P c (x,x) = 0 for all x € T", and hence h(y{t),y{t)) = 0. It 
follows from (2.11) that 7 (t) £ A , and thus (7(t),7(t)) € A, as we wanted to 
prove. □ 

We refer the interested reader to [21, 29] for more details on Mather-Mahe 
theory. 
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The minimal action and convex billiards 



In the classical spectral problem one would like to calculate or estimate the 
spectrum of a given Riemannian manifold. The inverse spectral problem asks 
the opposite question: how much information about the underlying manifold 
is encoded in its spectrum? In other words, to what extent is the geometry 
determined by the spectrum? Or, formulated by Kac [49] in his famous title: 
can one hear the shape of a drum? 

The term “spectrum” can have different meanings here. Either it stands for 
the eigenvalue spectrum of the Laplacian, or it means the length spectrum of 
the geodesic flow. These two interpretations are not completely independent, 
and there are subtle relations via the Poisson relation; see, for instance, [18, 
38]. It is known by now that the answer to Kac’ question is negative. There 
are whole families of isospectral, non-isometric manifolds; see, e.g., the survey 
in [35] . An in-depth survey on positive results concerning the inverse spectral 
problem has recently been given by Zelditclr in [105]. 

In the following, we want to consider the inverse spectral problem for 
strictly convex domains 17 in R 2 . The Laplace spectrum with Dirichlet bound- 
ary conditions describes the sound you hear when you beat the “drum” 17. The 
notion of length spectrum needs some explanation because 17 has a bound- 
ary. By definition, geodesics in a bounded domain 17 are geodesics (in our 
case, straight lines) that get reflected at the boundary according to the law 
“angle of reflection = angle of incidence” . Such geodesics are often called bro- 
ken geodesics. Then the length spectrum consists of the lengths of all closed 
(broken) geodesics, together with their multiples. 

One way to think of broken geodesics is to image a room 17 lined with 
mirrors on its side <912; broken geodesics are light rays in this mirrored cham- 
ber. Another way of looking at broken geodesics is to think of sound travelling 
inside the room 17. Finally, and this is the point of view that we will adopt, 
broken geodesics are the trajectories of a billiard “ball” (which is just a point 
rather) being played inside the billiard table 17 and going around without 
friction. 



K.F. Siburg: LNM 1844, pp. 37-57, 2004. 
© Springer- Verlag Berlin Heidelberg 2004 
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3 The minimal action and convex billiards 



The main question is whether we can recognize the domain 17 , respectively 
its boundary curve, from the knowledge of its length spectrum. One way to 
attack this problem is to construct length spectrum invariants (LS-invariants, 
for short) and to relate them to geometry. Difficulties may arise in a twofold 
way - to prove that a certain geometric quantity is an LS-invariant, or to give 
a geometric meaning to some known LS-invariant. 

The crucial observation is that, for planar convex domains, one can find 
coordinates such that the billiard ball map is a monotone twist map. Moreover, 
the length of a closed geodesic is, up to sign, the action of the corresponding 
orbit; hence length maximizing geodesics correspond to minimal orbits. This 
observation allows us to apply techniques from Aubry-Matlrer theory. We 
will see that the minimal action is invariant under continuous deformations of 
the domain that preserve the length spectrum. In addition, many geometric 
quantities — such as the lengths and Lazutkin parameters of convex caustics — 
can be read off from the minimal action. Note that, once we know that the 
minimal action is an LS-invariant, the proof that some geometric quantity 
obtained from it is also an invariant becomes trivial. Finally, we show that 
the asymptotics of the minimal action is determined by the Dirichlet spectrum 
of ft. 



3.1 Convex billiards 

Let 17 be a strictly convex domain in R 2 with C 3 - boundary df2. As for the 
regularity of <917 , we just point out the surprising observation by Halpern [40] 
that a C 2 -curve <912 may produce a geodesic flow which is not defined for 
all times. Our C 3 -condition guarantees the completeness of the flow. In this 
chapter, we will always assume that the length Z(<912) of the boundary curve 
is normalized to 1. 

A broken geodesic in 17 is completely determined by its reflection points, 
together with the angles of reflection. The map 

<j) : S 1 x (0, 7r) — > S 1 x (0,7r) 

(s 0 , V’o) I-* (si,if>i) 

that associates to a pair (s, ^>)=(arclength on <912, angle with the positive 
tangent) the corresponding data at the next reflection, is called the billiard 
map associated to 17. Let us denote by 

h(s,s') = -\P(s)-P(s')\ 

the negative Euclidean distance between two points on <917. Elementary ge- 
ometry (see Fig. 3.1) shows that 



f|(so,si) = cosi/jo 
0(s o ,si) = -cosVh 



(3.1) 
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Fig. 3.1. Proof of (3.1) 



as well as dsi/dipo > 0. 

Thus, if we lift everything to the universal cover and introduce new coor- 
dinates (x, y) = (s, — cos VO eRx (—1, 1), we have 

yi da;i - y 0 dx 0 = cl/i(a;o,;ri) 

as well as 

dxi _ 1 ds! > Q 

dy 0 sin 0o dip 0 

This proves the following proposition. 

Proposition 3.1.1. In the coordinates (x,y), the billiard map 
0:Rx (— 1, 1) — > K x (-1,1) 

is a twist map, with the negative Euclidean distance being a generating func- 
tion. 

Moreover, we may extend the billiard map to the closed strip S 1 x [—1, 1] 
by fixing the boundaries pointwise. 

A periodic orbit of the billiard map (f> corresponds to a closed (broken) 
geodesic inside the domain 17. In order to distinguish topologically different 
closed geodesics, we associate to each periodic orbit its rotation number as 
defined in Def. 1.2.2. Let us give a more geometric definition here. 

Definition 3.1.2. The rotation number of a periodic billiard trajectory (re- 
spectively, a closed broken geodesic) is the rational number 
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m winding number ^ h 
n number of reflections ’ 2 J ’ 

where the winding number m > 1 is defined as follows. Fix the positive orien- 
tation of dFl and pick any reflection point of the closed geodesic on dFl; then 
follow the trajectory and measure how many times it goes around <917 in the 
positive direction until it comes back to the starting point. 

Note that we restrict ourselves to rotation numbers less than or equal to 
1/2, since a closed geodesic with rotation number u> can be seen as one with 
rotation number 1 — u>, traversed in the backward direction; see Fig. 3.2. 



(fixed) orientation 




Fig. 3.2. Closed geodesics of rotation number 1/5 and 4/5 



It was G.D. Birklroff who introduced convex billiards as a conceptually 
simple, yet mathematically complicated, dynamical system. Applied to convex 
billiards, Birklroff’s theorem (Tlrm. 1.2.4) shows that for every to / n £ (0, 1/2] 
in lowest terms, there are at least two closed geodesics of rotation number 
to/ n. In fact, one of them is an inscribed n — gon with winding number to, that 
maximizes the perimeter amongst all such n-gons; the other one corresponds 
to a saddle point of the length functional. 

Example 3.1.3. Consider 2-periodic billiard trajectories in a strictly convex 
domain FI. Geometrically, they correspond exactly to diameters of 17, i.e., 
segments that meet the boundary at a right angle at both ends. Therefore, 
we have n = 2 and to = 1. Already the ellipse is an example with precisely 
two 2-periodic billard trajectories, showing that the lower bound in Birklroff’s 
theorem is sharp. 

Definition 3.1.4. The marked length spectrum of a strictly convex domain 
17 is the map 
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MC{Q) :Qn (0,^] 

that associates to any m/n in lowest terms the maximal length of closed 
geodesics having rotation number m/n. 

The length spectrum of i? is defined as the set 

£(17) := N {lengths of closed geodesics in 17} U N l(df2). 

Note that, due to Birkhoff’s theorem, the marked length spectrum is a 
well defined map. Moreover, l(dfi) = 1 by our standing assumption that the 
boundary length is normalized to 1. 

The length spectrum contains information about all closed geodesics, al- 
beit in an “unformatted” form. In contrast, the marked length spectrum does 
give the labelling by the rotation number but only for the closed geodesics of 
maximal length. 

Have you ever visited the great basilica in Rome and stood in its huge dome 
(42m in diameter)? If you are inside the domed roof and try to communicate 
with a friend on the other side of the dome. Rather than shouting into the air, 
get close to the circular wall and whisper along the wall - you will be heard 
clearly on the other side. 

This is the effect of what is usually called a “whispering gallery”. The 
sound waves get reflected and travel along the wall, always staying close to it. 
In the context of billiards, such a whispering gallery is called a caustic. 




Definition 3.1.5. Let f2 be a strictly convex bounded domain in R 2 . A convex 
caustic is a closed C 1 -curve in the interior of 17, bounding itself a strictly 
convex domain, with the property that each trajectory that is tangent to it 
stays tangent after each reflection; see Fig. 3.3. 
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More generally, a caustic is defined a a continuous curve inside 17 with the 
above property; it need not be differentiable, nor bound a convex domain. For 
our purposes, however, it is sufficient to restrict ourselves to convex caustics. 

Example 3.1.6. The simplest example is given by the disk 17 of perimeter 1. 
In the original coordinates (s, ip), the billiard map is given by 

(siiV’l) = ( s 0 + — > V’o)- 

7 r 

Since it preserves the angle, it is an integrable twist map in the coordinates 
(x,y) = ( s , — cosijj). Its phase space is foliated by invariant circles; the phase 
portrait is shown in Fig. 3.4. Moreover, the disk 17 is foliated by concentric 
circles, each of which is a convex caustic for the circular billiard; see also 
Fig. 3.4. 




(a) Caustics are concentric cir- 
cles 



y 




Fig. 3.4. The billiard inside a disk 



Example 3.1.7. The next simple example is the billiard inside an ellipse. It is 
known that this system possesses an integral, i.e., there is a (non-constant) 
quantity which is preserved along a trajectory [24] . This integral has the fol- 
lowing geometric interpretation. 

Consider the billiard inside an ellipse with foci i 7 ), Then each trajectory 
either 

1. always intersects the open segment between the two foci, or 

2. always passes through the two foci alternately, or 

3. never intersects the closed segment between the foci. 

In fact, each trajectory, which does not pass through a focal point, is always 
tangent to precisely one confocal conic section, either a confocal ellipse (in 
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which case the trajectory never intersects the segment between the foci) or 
the two branches of a confocal hyperbola (where the trajectory always in- 
tersects the segment between the foci); see Fig. 3.5. The eccentricity of the 
corresponding conic section, for example, can be taken as an integral for the 
elliptical billiard. For proofs and further remarks, the reader may consult 
[24, 101]. 

Thus, the confocal ellipses inside an elliptical billiard are convex caustics 
in accordance with Def. 3.1.5, so the elliptical billiard is foliated by convex 
caustics (up to the segment between the foci). The branches of the confocal 
hyperbolae can then be seen as caustics in the more general sense mentioned 
above. 




(a) Caustics are confocal ellipses and hy- 
perbolae 



y 




Fig. 3.5. The billiard inside an ellipse 



The phase portrait of an elliptical billiard is also shown in Fig. 3.5. Al- 
though it looks like the phase portrait of the pendulum (Fig. 1.3), the dynam- 
ics are quite different. The points (0, 0) and (1/2, 0) and its translates do not 
represent equilibrium points anymore, but belong to the two-periodic orbits 
corresponding one of the half-axes of the ellipse, and similarly for the other 
half-axis. Their rotation number is 1/2, which implies that the islands are not 
fixed, but “wander”: they are mapped onto each other. 

Bounding the islands we see separatrices, corresponding to the orbits 
through the foci. The invariant curves above and below the seperatrices rep- 
resent the orbits not intersecting the segment between the foci (i.e., being 
tangent to confocal hyperbola). 

As an aside, we mention here that a famous conjecture, usually attributed 
to Birklroff, states that the elliptical is, in fact, the only convex billiard with 
an integral. 
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Let us return to the general case of a convex billiard 17. Suppose for a 
moment that the billiard possesses a convex caustic c. Then one can associate 
the following two parameters to c : 

1. its rotation number w £ (0,1/2), defined as the rotation number of the 
circle homeomorphism on c induced by the geodesic flow via the points of 
tangency; 

2. its length 1(c). 

It turns out that there is a third parameter associated to a convex caustic, 
the so-called Lazutkin parameter. 

Definition 3.1.8. Let ft be a convex billiard with a convex caustic c. Then 
the Lazutkin parameter of c is defined as 

Q(c) = \A-P\ + \P-B\-\AB\, 

where P is any point on dfi and A, B £ c are the points of tangency of c seen 

from P; see Fig. 3.6. Moreover, \AB\ denotes the length of the caustic’s part 
from A to B, where we have oriented the caustic according to the geodesics 
touching it. 




Fig. 3.6. The Lazutkin parameter of a convex caustic 



In fact, if c is not a caustic but just any closed convex curve inside 17, the 
Lazutkin parameter can be defined in the same manner but may depend on 
the point P £ dfi. It is independent of P if, and only if, c is a caustic [55, 1]. 
Therefore, the Lazutkin parameter of a caustic is well defined. 

What is the relation between (convex) caustics of a convex billiard 17 and 
invariant circles for the corresponding billiard map <f)l Certainly, to a convex 
caustic in 17 corresponds an invariant circle for the billiard map, i.e. a simply 
closed, homotopically nontrivial curve T in S 1 x (—1,1) with <j>(r) = r. The 
converse, however, is not entirely true. By a theorem of Birklroff (see [94] 
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and the references therein), invariant circles of twist maps are graphs and 
therefore do give rise to caustics; but these caustics need neither be convex 
nor differentiable. 

Finally, one may ask whether convex caustics exist for some arbitrarily 
given convex billiard. Lazutkin proved [55] that the billiard map of a convex 
billiard with sufficiently smooth boundary possesses a Cantor set of invariant 
circles near the boundary, a result which is based on Moser’s twist theorem 
on invariant curves [75] (hence the condition on being sufficiently smooth). 
Fortunately, invariant circles near the boundary always correspond to caus- 
tics which are convex C-curves, so every convex billiard with C°°- boundary 
possesses (uncountably many) convex caustics according to Def.3.1.5. 



3.2 Length spectrum invariants 

In this section, we start to investigate which geometric data of a strictly 
convex domain IlcI 2 are determined by its length spectrum. Recall from 
Def. 3.1.4 that the length spectrum is the set C(f2) consisting of all multiples 
of lengths of closed (broken) geodesics in 17. The marked length, on the other 
hand, is the map associating to every rational m/n the maximal length of a 
closed geodesic of rotation number m/n. In order to study the length spectrum 
one looks for length spectrum invariants, i.e. , quantities that depend only on 
the length spectrum; see Def. 3.2.1 below. In particular, we are interested in 
invariants that carry some geometric information about the domain 17. 

In general, it is not clear to what extent the marked length spectrum is 
determined by the length spectrum, or vice versa. For this reason, we make 
the following definition. 

Definition 3.2.1. A quantity (number, function, etc.) is called a (marked) 
length spectrum invariant if for any two strictly convex domains 17 0 , L?! C M 2 
with the same (marked) length spectrum this quantity is the same. 

We will write LS-invariant for length spectrum invariant, and MLS- 
invariant. for marked length spectrum invariant. 

As mentioned before, the notions of LS-invariant and MLS-invariant may 
differ. For continuous deformations of smooth (C°°) domains, however, one 
has the following result. 

Proposition 3.2.2. Suppose 17 s ,s € [0,1], is a continuous family of strictly 
convex domains with C°° -boundaries such that £(17 S ) = C(f2o) f or o,ll s. Then 

MC{n s ) = MC({2 0 ) 



for all s. 
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Proof. Pick a rational rotation number m/n and a closed geodesic in 12 s of ro- 
tation number m/n having maximal length. We claim that the corresponding 
value MC(f2 s )(m/n) is independent of s. 

The chosen geodesic corresponds to a periodic orbit of the corresponding 
twist map, respectively, to an (n + l)-tuple with £o < £1 < • • • < = Co + m 

which is a minimum of the finite action functional 

n 

H s (foi • ■ • > £n+l) = ^ ' ^s(£i> £i+l)> 

»= 0 

where the generating function h s is nothing but the negative length function 
in fi s . By Sard’s Theorem, the critical values of H s , i.e., the values of C(I2 S ), 
form a set of Lebesgue measure zero. In general, the maxima of a family of 
smooth functions depend continuously on the parameter. In our case, they lie 
in a set of Lebesgue measure zero, so they must stay constant. 

This means that the maximal length of closed geodesics having a fixed 
rotation number m/n does not depend on s. Hence the functions A 4C(I2 S ), s € 
[0,1], are all the same. □ 

Corollary 3.2.3. For continuous deformations of smooth strictly convex do- 
mains I2 S , MLS-invariants are also LS -invariants. In other words, if a certain 
quantity stays invariant under deformations with A 4C(f2 s ) = Ai£(I2o), it is 
also invariant under deformations with C(I2 S ) = C(I2q). 

We have seen in the previous section that the billiard map associated to a 
convex domain 12 is a monotone twist map on S 1 x (—1, 1) generated by h, the 
negative distance between points on df2. Therefore, closed geodesics of maxi- 
mal length correspond to minimal orbits, and the marked length spectrum is 
essentially nothing but the minimal action; more precisely, 

o(— ) = -~MC{Q)(—) (3.2) 

n n n 

for every m/n € (0, 1/2] in lowest terms. 

This simple observation turns out quite fruitful because it implies the 
following principle. 



Main Principle. Every quantity that can be calculated from the minimal 
action is, by tautology, a marked length spectrum invariant. 



Remark 3.2./. One of the main advantages of this principle is that the actual 
proof that a certain quantity is a MLS-invariant becomes trivial, once it is 
clear how the quantity can e calculated from the minimal action; we will see 
applications of this remark in the following sections. 

Vice versa, we see from the identity (3.2) that every MLS-invariant must 
be hidden in the minimal action — the only question is how. Let us formulate 
the above principle as a theorem. 
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Theorem 3.2.5. The minimal action a : [0, 1] — > R. for the billiard map 
of a strictly convex C 3 -domain Q C M 2 is a complete MLS-invariant, i.e., 
A4£(]?o) = A4£(J?i) if and only if ao = oi- 

Moreover, a is a strictly convex function on [0, 1], symmetric with respect 
to the point 1/2, and three times differentiable at the boundary points with 
a'(O) = -l{df2) = -1. 

Proof. The first assertion follows from (3.2) and the continuity of a. The strict 
convexity is contained in Proposition 1.3.7; the symmetry property is obvious. 

For the last part, we make use of a special choice of coordinates near the 
boundary {ip = 0} which is due to Lazutkin [55, 56]. Namely, let us introduce 

£ = C J p~ 2/3 (r) dr , 77 = 4Cp 1 / 3 (s) sin ^ 

where p € C 1 (S 1 ,K) is the radius of curvature and C = ( f 4 p~ 2 / 3 )~ 1 . Then 
{ip = 0 } corresponds to {77 = 0 }, and the invariant symplectic form sin ip dip A 
ds takes the form 

^ 7/ cfy A rM) . 

Moreover, one calculates that 

= (Z + v,v) + o(v 2 ) 

as 77 —7 0. Hence (p is a perturbation of the integrable twist map (£, 77 ) 1 — > 
(/ + 77 , 77 ), and we obtain for the generating function 1 

M£o, £ 1 ) = (£1 - £o) 3 + C>((£i - £o) 4 ) 

as £ 1— £0 — ► 0. Prop. 1.3.8 implies that the minimal action — in the coordinates 
£, 77 — can be written as 

a Lv(u) = + ^(w 4 ) 

as to — > 0 . 

Transforming back to the coordinates x, y means adding the linear term 
— U) to the action [67] because the cohomology class [ydx — ( 8 C 3 ) - 1 ? 7 2 df] is 
— 1 ; the latter is due to the fact that the boundary {ip = 0 } corresponds to 
{y = —1} respectively {77 = 0}. Summarizing, we have for the minimal action 

«(u) = -w+— w 3 + d(w 4 ) (3.3) 

as u> — > 0. By symmetry, an analogous formula holds near u> = 1. □ 



1 Note that the exponent in the first term is 3, due to the nonstandard form of the 
symplectic form. 
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Recall from (1.7) that the convex conjugate of a is defined as a* (I) = 
max u \u>I — a(u>)\ with (a*)' (a' (lu)) = u>. Hence, in view of Theorem 3.2.5, 
the domain of definition of a* is the interval [—1, 1]. By way of illustration, 
let us calculate the minimal action and its convex conjugate in the simplest 
example. 

Example 3.2.6. Take fi to be the disk of perimeter 1 in R 2 . The billiard map 
is integrable and given by (si,ipi) = (so +V , o/ 7 b^o) with generating function 

h(s, s') = sin7r(s — s'). 

ir 

The whole phase space is foliated by invariant circles {y = const.}, and 
the minimal action is just 

a(u>) = — — sin7rw. 

7 r 

In view of the identity sin arccos(— x) = \/l — x 2 , its convex conjugate is 

a* (I) = — ^ arccos (— /) • I + s/l — I 2 ^ € [0, 1] 

for / £ [—1, 1], The graphs of a and its convex conjugate a* are depicted in 
Fig. 3.7. 




Fig. 3.7. The minimal action for the disk and its convex conjugate 

The asymptotics for uj — > 0 and I — 1, respectively, are as follows: 

a(uj) = -u+ yw 3 - + 0(uj 7 ) 

«*(/) = ^ (| (i + 1) 3/2 + + 1) 5/2 + Jo (/ ■ + 1) 7/2 ) + 0((I + l) 9 / 2 ) 

a can be extended to an odd smooth function on R. a* has a singularity of 
order 3/2 at I = —1; the function (ct*) 2 / 3 is smooth. 
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If we ask whether the minimal action can be recovered from the length 
spectrum instead of the marked length spectrum, we have the following result 
for continuous deformations. 

Corollary 3.2.7. For continuous deformations of smooth strictly convex do- 
mains, a is an LS-invariant function on [0, 1] . 

Proof. Immediate from Corollary 3.2.3. □ 

Finally, we point out that the minimal action is actually invariant under 
arbitrary symplectic (i.e. , area-preserving) coordinate changes in the phase 
space. For, by (3.11), the actions of periodic orbits are symplectically invari- 
ant, and they determine a. 



3.2.1 Classical invariants 

What kind of geometric data can be recovered from the (marked) length 
spectrum? It is well known, for instance, that the perimeter and the diameter 
of the domain fl are invaraints of the spectrum, a fact which follows readily 
from the Main Principle. 

Proposition 3.2.8. The boundary length l{dfi), the diameter diamfi, and 
the curvature integral f p~ 2 / 3 are MLS -invariants. 

Proof. This is an immediate consequence of our Main Principle, together with 
(3. 3), since l(df2) = —a'(0), diaml? = — a(l/2), and (/ p -2 / 3 ) 3 = 4a"' (0). □ 

Remark 3.2.9. In view of Cor. 3.2.7, the above quantities are LS-invariants 
for continuous deformations of smooth strictly convex domains. 

Are there other geometric data hidden in the minimal action? The follow- 
ing theorem shows that the parameters connected to a convex caustic can be 
read off from the minimal action. 

Theorem 3.2.10. Let Q be a strictly convex C 3 -domain in R 2 , and suppose 
that c w is a convex caustic of rotation number to. Then the length of c w and 
its Lazutkin parameter are given by 

l( c u ) = -a'(ur) 



and 

Q(c u ) = a*(a\u)). 

Proof. Call T(s) and N(s) the unit tangent vector and unit inward normal, re- 
spectively, at a point P(s) € dfi, and set U(s) = cos rp(s)T(s) +smrp(s)N(s); 
here, ip(s) £ (0, 7t/2) is the unique angle such that the ray from P(s) having 
direction U(s) touches c w . Then there is a function t(s) such that 
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A(s) = P(s ) + t(s)U(s) G 



and 

A(s) =T(s) + t(s)U(s) + t(s)U (s) || U (s). 

Since A(s) 0, we can write 

|i(s)| = (A(s),U(s)) = COS lf(s) + f(s) 

so that 

T(c u ) = / |A| ds = / cos ipds = — / ydx——a'{ijj) 

J s 1 J s 1 Jr, 

by Proposition 1.3.7. 

By definition of the Lazutkin parameter (see Fig. 3.6), we have 



N 



N+l 



(N + l)Q(c w ) — |Ai — Pi| + |P — P+i| + \Pn+i ~ Pjv+i| — | AiBi 

i- 1 i=l 

for TV > 1. Hence 



1 N 1 N ^ 
^ N^ lPi ~ Pi+l1 " a!™, TV £ ^ 

2 — 1 i— 1 

= —a(cj) - u>L( c^) 

= wa'(w) — a(a>) 

= a*(a'M) 



□ 

Corollary 3.2.11. 77ie length and Lazutkin parameter of a convex caustic are 
MLS-invariants, respectively, LS-invariants under continuous deformations. 

Remark 3.2.12. The result that lengths and Lazutkin parameters of convex 
caustics are spectrally determined is due to Amiran [2]. The above proof, 
however, is new. 

Remark 3.2.13. In view of Theorem 3.2.10, one might call — a ' and a* the 
generalized “length” and “Lazutkin parameter”, even if there is no convex 
caustic of the corresponding rotation number. 

Given a convex caustic c, one can reconstruct <917 by wrapping a string of 
length 1(c) + Q( c) around c, pulling it tight, and going along c. These “string 
length parameters” of convex caustics are, of course, also MLS-invariants of 
the domain 17. 

The following result shows how analytical properties of a translate into 
geometric properties of 17. 
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Proposition 3.2.14. Let fl C R 2 be a strictly convex C 3 -domain with as- 
sociated minimal action a : [0,1] — > R. Then the following statements are 
equivalent: 

1. a is differentiable at the point 1/2; 

2. dfl is a curve of constant width. 

Proof. The differentiability of ct at 1/2 is, in view of Thm. 1.3.7, equivalent to 
the existence of an invariant circle consisting entirely of minimal 2-periodic 
orbits, i.e., diameters of the same length. □ 

Another well known fact is the spectral rigidity of the disk, i.e., that the 
disk is completely determined by its spectrum. In our context, we have the 
following rigidity result. 

Theorem 3.2.15. Let 17 C R 2 be a strictly convex C 3 -domain with associated 
minimal action a : [0, 1] — > R. Then the following statements are equivalent: 

1. a is differentiable; 

2. a is analytic; 

3. fl is a disk. 

Proof. Since 3. => 2. => 1. in view of Example 3.2.6, it suffices to show 1. => 3. 

According to Thm. 1.3.7, the differentiability of a at the point 1/2 is 
equivalent to the existence of an invariant circle consisting of minimal periodic 
orbits. Taking limits of these curves, we obtain invariant circles for all rotation 
numbers, consisting entirely of minimal orbits. 

We claim that they foliate the phase space. Indeed, if there was a gap, its 
boundary curves would necessarily have the same rotation number (otherwise, 
there would be rotation numbers without invariant circles). But this is impos- 
sible, due to the graph property of the set of minimal orbits in Theorem 1.3.4. 

Now the assertion follows from Bialy’s result [9] that the only billiard 
whose phase space is foliated by invariant circles is a circular one. □ 

Remark 3.2.16. We conjecture that the statement “a is differentiable near the 
point 1/2 => fl is a disk” is true, which would mean that the differentiability 
of a minimal action, stemming from a strictly convex billard, near 1/2 already 
implies its global differentiability. 

Corollary 3.2.17. Suppose fl is a strictly convex C 3 -domain with the same 
marked length spectrum as the disk of the same perimeter. Then fl is a disk. 

Remark 3.2.18. In view of Corollary 3.2.3, we have the following rigidity result 
for deformations: any continuous deformation of the disk (inside the class of 
smooth convex domains) that preserves the length spectrum must be trivial. 

Remark 3.2.19. Corollary 3.2.17 has an analogue in differential geometry: a 
Riemannian 2-torus, having the same marked length spectrum as a flat torus, 
is flat [7, Thm. 6.1]. 
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3.2.2 The Marvizi— Melrose invariants 

So far, we have considered the most general case of domains with C 3 - 
boundaries where we do not know, for instance, whether there are any con- 
vex caustics at all. If we assume smooth enough boundaries, however, then 
Moser’s twist theorem [75] can be applied and guarantees the existence of a 
Cantor family of convex caustics that accumulate at the boundary dfi; this 
was first proven by Lazutkin [55] . From now on, we assume that 17 is smooth 
(C°°). Then there is the following KAM-theorem by Kovachev and Popov 
[52, Thm. 2] which is based on the work of Poschel [89]. 

Theorem 3.2.20. Let Li be a smooth strictly convex domain of unit boundary 
length. Then there are symplectic coordinates (0,1) near {y = —1} <-» {I = 
— 1} such that the billiard map 4> '■ (Oo,Iq) > (0\,li) is generated by 

5(0o, h) = 0 o h + K(hf/ 2 + R(0 o , h), 



i.e., 

f/o = d\ S = I\ + d\R 

\0i = d 2 S = 0 o + |i^(/ 1 ) 1 / 2 AT , (/ 1 ) + d 2 R 

Here , K G C°°(R,M.) with K(- 1) = 0,iC(-l) > 0, and R G C^R^R) 
is 1 -periodic in the first variable. Moreover, there exists a Cantor set C* C 
[-1,-1 + e*) with —1 G C* , where e* > 0 is some small number, such that 
R = 0 on 1 x C*. 

We see that the perturbation term R vanishes on IxC* with all its 
derivatives. Each curve R x {I}, I G C* , gives rise to an invariant circle for 
the billiard map on which it is conjugated to the rigid rotation 

(0, 1) ^ (0 + oj, I) (3.5) 

with oj = 3/2K(I) 1 / 2 K i (I). Since these invariant circles lie near the boundary, 
they correspond to a Cantor set of convex caustics near and accumulating at 

on. 

In view of (3.5), 0 o Ii + a*(I\) generates (f : (0 o,Iq) > (0\,I\) on the 

Cantor set RxC* of invariant circles accumulating at R x {— 1}. We also know 
from Theorem 3.2.20 that </> is generated by S = 0 o I 1 + A'(/ 1 ) 3 / 2 + R(6 0 , Ii), 
where K is smooth and R vanishes with all its derivatives on RxC*. Since 
a*(— 1) = 0 = K(— 1), we must have that 

a*(I) = K(I ) 3 /2 (3.6) 

for all / G C*. Thus, a* yields a generating function for (f>, restricted to a 
Cantor set of invariant KAM-circles. 

The smooth function K can be written as 



K(I) = K\- 1) •(/+!) + 0((I + l) 2 ) 



(3.7) 
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with K'(— 1) > 0. Combining (3.6) and (3.7) we obtain that 
a (tv) = aiu) + a^co 3 + O(to 5 ) 
is smooth on the Cantor set 

C= ^K{lf/ 2 K\I) | JeC’} 



(3.8) 



containing 0. 

Let us return to the original question which geometric data can be re- 
covered from the (marked) length spectrum. Concerning the lengths and 
Lazutkin parameters of convex caustics — whose existence is now guaranteed 
by Thm. 3.2.20 — Cor. 3.2.11 immediately implies the following assertion. 

Theorem 3.2.21. The function tv i— > (l(c u ), Q(c u )), defined on a Cantor set 
in [0, 1] containing 0 and 1, is an MLS-invariant, respectively an LS-invariant 
under continuous deformations. 

Remark 3.2.22. Again, we point out that Thm. 3.2.21 is not new: it is the 
main result in [88]. Popov’s proof, however, relies on the existence of invariant 
KAM-circles, and a good part of it is hard analysis. Our approach is simpler 
and shows that the invariance property of the caustic parameters is definitely 
not a phenomenon inside KAM-theory. 

In the smooth case, Theorem 3.2.10 also implies a functional dependence 
of 1(c) and Q(c). 

Theorem 3.2.23. There is a formal power series expansion 

1= 1 + bkQ 2k/3 

k> 1 

as Q — > 0, whose coefficients are MLS-invariants, respectively LS-invariants 
under continuous deformations. 

Proof. l( c u ) = —a'(tv) and Q( c w ) 2 / 3 = a* (a' (tv )) 2 / 3 = K(a'(tv )) are smooth 
functions on C with K'(— 1) > 0. The claim follows from the implicit function 
theorem. □ 

Remark 3.2.24- Cor. 3.2.23 was first stated explicitly in [2, (3.1)]. It follows, 
however, already from formulae obtained by Lazutkin, namely (1.11) and 
(1.12) in [55], where he expresses the length and his parameter in terms of the 
rotation number. 

Let us go on and see if there are other spectral invariants hidden in the 
minimal action. In 1982, Marvizi and Melrose [63] defined a sequence of MLS- 
invariants by investigating an integrable approximation of the billiard map 
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near the boundary. The Marvizi-Melrose invariants are given by the asymp- 
totics of a so-called interpolating Hamiltonian. By definition, an interpolating 
Hamiltonian is a smooth function ( on S 1 x [—1, 1], whose time-^^-map is 
the billiard map, up to a diffeomorphism that fixes the boundary to infinite 
order. The integral invariants are then defined as the Taylor coefficients of 



J(r) 



1 

C'(C -Hr)) 



(3.9) 



at r = 0 for any interpolating Hamiltonian £ in action-angle-variables. 

We claim that, in the context of minimal action, we can find a quite simple 
interpolating Hamiltonian. Indeed, from (3.4) and (3.6), we see that 






2/3 



is an interpolating Hamiltonian, at least on a Cantor set containing the bound- 
ary. Since the Taylor coefficients of (3.9) at the boundary point r = 0 only 
depend on the behaviour of <£(/) on any sequence accumulating at 0, they only 
depend on the Taylor coefficients of a* at the corresponding boundary point 
— 1. Summarizing, we obtain the following result. 



Theorem 3.2.25. The integral invariants of Marvizi and Melrose are alge- 
braically equivalent to the Taylor coefficients of a* at —1. 

Marvizi and Melrose [63] go further and prove an asymptotic formula for 
the lengths of closed geodesics which is then shown to be spectrally deter- 
mined. In order to integrate this formula into our context, we consider a 
closed geodesic g mn of rotation number m/ n. Its length is given by 



l{.9mn) — ^ ' h. 

Qmn 



(3.10) 



We want to rewrite this in ( 6 , J)-coordinates and relate it to the generating 
function S(9q,Ii) from Theorem 3.2.20. We have the following transforma- 
tions: 



d> : ( 9 , 1) i— > (a;, y) with T>* ( y dx) — I d9 = dH 
q/ — o f)oT> 

9\ dl\ T Iq d9o = dS with S = 9qIi 4- K[Ii)^^ ~ -T R{9o, I \ ) 

A straightforward calculation shows that the generating function transforms 
according to the formula 

T'*(Id9)-Id9 = d(ho$ + H-HoT>). (3.11) 

On the other hand, we can write the left hand side as 
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I\ dd\ — Iq dOo — ( — 61 dli — Io dd 0 ) T (0i dl\ I\ dO i) 

= - ds+ (f 

= dS* 

From this we conclude that 



— h o <P = H — Hod'— S* + const. 

Summed over a closed orbit, the term H — H o A adds to zero, so (3.10) yields 



-l(9mn) = - V-S* 
n n z ' 

Qmn 



const. 



n 



Since S* = (A'(/) 3 / 2 )* + I dR/dl — R has the same Taylor series at the 
boundary point as (AT (A) 3 / 2 )*, which coincides on S 1 x C with (a*)* = a, 
we see that l(g mn ) / n has the same Taylor series for n — > oo as —a. Loosely 
speaking, all closed geodesics become minimal as they approach <912. More 
precisely, 

lidmTi) = TYl T~ 'y ) C-mkW’ 
k> 1 

where 

rn 2k-\-l 

Cmk = ~ (2fcTT)I a{2k+1) ^- ^ 

This formula implies the following 



Theorem 3.2.26. If ft is smooth, the asymptotics of the length spectrum as 
n — > oo are equivalent to that of the minimal action as lo — > 0. 



3.2.3 The Gutkin Katok width 

Another application of Theorem 3.2.5 concerns regions in ft which are free of 
convex caustics. Gutkin and Katok gave estimates for their area in terms of 
the geometry of 12. In particular, they proved [39, Prop. 1.3] that a convex 
caustic c u with rotation number lo cannot lie too far from the boundary: 

max d(P, dI2) < yMiam.!? • Q(c u ). 

Pe 

See Fig. 3.8 for an illustration. 

We will see that the Gutkin-Katok width i/diaml2 • Q( c w ) has an MLS- 
invariant interpretation. Indeed, in view of Prop. 3.2.8 and Theorem 3.2.10, 
this number can also be written as 

\J diaml? • Q(c w ) = \/-a(l/2) ■ a*{a'{u)). 

According to our Main Principle, this proves the following result. 
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Fig. 3.8. The Gutkin-Katok width 



Theorem 3.2.27. Suppose 12o,l7i are two strictly convex C 3 -domains with 
the same marked length spectrum. Then, for a fixed rotation number to, every 
convex caustic in l7o, respectively fl\, is contained in a strip around <9 17 o, 
respectively dfl\, of one and the same width. 
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As already mentioned in the introduction, there is a relation between the 
length spectrum £(17) of a strictly convex domain and the spectrum of the 
Laplacian with Dirichlet boundary conditions: 



A u = A 2 u in 17 
u = 0 on <917 



(3.13) 



Indeed, due to a Poisson relation for compact Riemannian manifolds with 
geometrically convex boundary [38] , the expression 

= Z cos At 

A 2 6specA 

is well defined as a distribution, which is smooth away from £(17) [3]. More 
precisely, if T > 0 is in the singular support of er then T lies in the length 
spectrum of 17. 

Conversely, whether some given T G £(17) belongs to the singular support 
of <t, depends on possible cancellations of singularities stemming from different 
closed geodesics of the same length. It is known [38] that T lies in the singular 
support of sigma if there is exactly one closed geodesic of length T, whose 
Poincare map does not have an eigenvalue 1. Marvizi and Melrose [63] showed 
that a much weaker non-coincidence condition on 17 suffices to conclude that 
almost all maximal lengths of geodesics having rotation number 1/n lie in 
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the singular support of a. Popov [88] generalized this to geodesics of rotation 
number m/n with m > 1, provided (m, n) is “near” the Cantor set described in 
(3.8). In particular, that non-coincidence condition is satisfied by all curvature 
functions in a C-neighbourhood of the constants. 

Thus, the values a(l/n), respectively a(m/n), with sufficiently large n are 
spectral invariants of the domain. This is also true for the coefficients cu- in 
(3.12), and hence for the Taylor coefficients of —a at 0. Therefore we can state 
the following result which is an analogue of [63, Thm. 7.4], 

Theorem 3.3.1. Suppose fi C R 2 is a smooth strictly convex domain with 
unit boundary length, such that 1 is not a limit point of lengths of closed 
geodesics having fixed rotation number m/n with m > 1. Then the Taylor 
series of the minimal action at 0 is completely determined by the Dirichlet 
spectrum (3.13). 

Applying Popov’s more general result, one can show that also the values 
of a on the Cantor set C (and hence all the caustic parameters L{ c) and Q{ c)) 
are spectral invariants under the non-coincidence condition in [88, (6.1)]. 

For an exhaustive treatise on the inverse (Laplace) spectral problem prob- 
lem we refer to the survey by Zelditclr [105]. 
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The minimal action near fixed points and 
invariant tori 



When investigating a dynamical system it is often a first important step to 
transform it into a form as simple as possible. If the system belongs to some re- 
stricted class these normal forms should be invariant under the corresponding 
invariance transformations. 

In classical mechanics, the class under consideration consists of Hamilto- 
nian systems, respectively symplectic transformations. We consider a symplec- 
tic map near a fixed point (in two dimensions), respectively, near an invariant 
torus (in arbitrary dimensions). Under certain nondegeneracy conditions on 
the linearization of the map , G.D. Birklroff constructed a normal form that 
is invariant under symplectic coordinate changes. This Birklroff normal form 
describes an integrable approximation of the original map. Its asymptotics at 
the fixed point define a set of symplectic invariants; for obvious reasons, we 
call these Birklroff invariants asymptotic. 

The goal of this chapter is the construction of a new local symplectic invari- 
ant that, in particular, includes the asymptotic Birklroff invariants. To do so, 
we will associate to the germ of a symplectic map at a fixed point /invariant 
torus the germ of its minimal action. This is symplectically invariant and, 
in contrast to the Birklroff normal form, reflects part of the dynamical be- 
haviour in a neighbourhood of the fixed point/invariant torus. We show that 
the classical Birklroff invariants are encoded in the minimal action as the Tay- 
lor coefficients of its convex conjugate. Moreover, in the integrable case, the 
minimal action determines the map completely, a fact which is also not true 
for the Birklroff normal form (unless the map is analytic). 

Symplectic mappings near a fixed point appear as Poincare section maps of 
a closed trajectory. We explain this for the geodesic flow on a compact surface. 
It turns out that the minimal action depends only on the length spectrum of 
the Riemannian manifold. Therefore, the minimal action is a new local length 
spectrum invariant for compact two-dimensional manifolds. 

In the final section we investigate the role of the minimal action near an 
invariant torus of a symplectic map. 
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4.1 The minimal action near plane elliptic fixed points 

We consider an area-preserving diffeomorplrism cj> of the plane in the vicinity 
of a fixed point. We can shift the fixed point into the origin, so we assume that 
<^(0) = 0. Moreover, we are only interested in the local behaviour of <f> near 
0. Therefore, we call two diffeomorphisms equivalent if they coincide on some 
open neighbourhood of the origin. In other words, we consider the equivalence 
class of (j>, i.e., we consider the germ of a symplectic diffeomorphism at the 
fixed point 0 £ (R 2 , 17) where 17 is some area form on R 2 . In the following, 
whenever we pick a representative </> : /7 — > R, we assume that U is a simply 
connected neighbourhood of 0; this is no loss of generality. Then, by Poincare’s 
Lemma, the symplectic form 17 is exact, i.e., there is a 1-form A with 17 = dX. 

Suppose for a moment that there is a point p £ U \ {0} whose iterates 
Pi = <t> l {p) exist for all i £ Z. We want to define the average action and the 
rotation number of the orbit (pi)iez- Since (f> is symplectic, the 1-form <f>* A— A 
on U is closed, hence exact: 



0*A-A = dS. 

Definition 4.1.1. Let <f> : U R 2 be a symplectic diffeomorphism, and A be 
a (local) primitive of the area form 17. 

Any function S : U —> R satisfying <j>* A — A = dS is called a generating 
function for </>. Any two generating functions differ by some additive constant, 
and we normalize S by setting 5(0) = 0; this makes the generating function 
unique. 

If (pi)i£ 7 , is an orbit of f> then the average action of (pf) is defined as 

A(( Pi )) := £ S( Pi ) 

i=-N 



if this limits exists. 

Remark f.1.2. As an aside, we remark that the normalization of the generating 
function is in accordance with setting H (t, 0) = 0 when H is a Hamiltonian 
whose flow generates <f> and leaves 0 fixed; for, then the generating function 
/ A — H dt vanishes at 0. 

Lemma 4.1.3. The average action does not depend on the choice of the 1 
form A. Moreover, it is invariant under local symplectic coordinate changes 
fixing the origin. 

Proof. We claim that the definition of the average action does not depend on 
the choice of the 1-form A. Indeed, taking another 1-form A' with d\' = dX = 
17, the closed 1-form A' — A is exact: A' — A = dF. The new (normalized) 
generating function is given by S' = S + (f*F — F. But the average of <f>*F 
over an orbit is the same as that of F, so the summation of <j>*F — F over the 
orbit vanishes. This proves our first claim. 
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Moreover, we claim that the average action is invariant under local sym- 
plectic coordinate changes <P fixing the origin. Such a transformation F is 
generated by some function F, due to Poincare’s Lemma. Then a short 
calculation shows that the generating function for F o <f> o F~ x is given by 
F*S + F — (F o (j>o F~ 1 )*F. Now the claim follows by the same argument as 
above. □ 

Next, we want to define the rotation number of an orbit (pi)jgz in U\ {0}. 
Roughly speaking, this is its average winding number around the origin. More 
precisely, we introduce polar coordinates on R 2 \ {0} = S 1 x (0,oo), and lift 
<f> : [/ \ {0} — > R 2 \ {0} to a map </> which is then defined on some strip 
in the universal cover R x (0,oo) of S 1 x (0, oo). Since <f> is an orientation- 
preserving diffeomorphism, <f> is a diffeomorphism of degree 1. Given an orbit 
( Pi ) of </> projecting onto (p,;) and a natural number N, we choose a curve 
r : [— N, iV] — > R x (0, oo) with r(i) = pi. Call r : [— N, N] — > R 2 \ {0} the 
projection of T, and close it up to a closed curve Gv by adding a “short” 
piece (whose lift upstairs lies inside one fundamental domain); see Fig. 4.1. 




Fig. 4.1. The rotation number of an orbit (pi)*ez 



Then we define the rotation number of (p,;) to be 

■= Jim [Fjv] G iLi(R 2 \ {0},R) 

N—>oo 

if this limit exists. Clearly, if we fix the lift q J, the class p does not depend 
on the particular choice of (pj) and Cv- Moreover, choosing a different lift (f> 
means adding the class [<j> — <j>] of the deck transformation <f> — cj). Finally, we 
have a canonical identification iii(R 2 \ {0}, R) = R if we take as generator of 
H i(R 2 \ {0},R) the class represented by the positively oriented unit circle. 




